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20.  Abstract  (continued 

295^K  and  the  vapor  pressure  are  calculated  throughout  the  col- 
lapse of  a supersaturated  metastable  vapor. 

An  effective  molecular  pair  interaction  energy  and  nearest 
neighbor  coordination  number  is  found  based  on  the  classical 
zero-point  enthalpy  of  sublimation  and  by  iteration  to  satisfy 
initial  equilibrium  limiting  constraints  at  273°k.  The  internal 
energy  redistribution  frequency  results  from  close  matching  of  the] 
experimental  pressure  distribution.  A method  of  determining  the 
cluster  size  at  which  macroscopic  properties  become  relevant  is 
described;  the  size  is  found  to  be  a cluster  of  cdx>ut  122  mole- 
cules at  273*K  based  on  a hard  sphere  model  collision  cross 
section.  For  the  water  vapor  dimer  an  equilibrium  constant^of 
4 . 0xl0"21(ra-3)  is  found  at  273*K  compared  to  3rlxl0“^l  (cm”"’) 
obtained"’*'  from  Keyes'('3)  data.  Two  models  for  water  vapor 
clusters  result  (from  satisfying  all  constraints  imposed  on  the 
theory)  which  are  given  in  terms  of  number  of  neeurest  neighbors, 
pair  interaction  energy  per  molecule  (2.9kcal/  . at  273*K>  and 

possible  structure  on  a cluster  by  cluster  basiS'^'for  a classical 
hard  sphere  model  and  for  a model  approximating  a Pauling  type 
clathrate. 
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I.  ABSTRACT 

The  kinetic  theory  of  cluster  formation  in  a condensing 
qas  proposed  by  Buckle attempts  to  calculate  the  course  of 
homogeneous  condensation  from  molecular  rather  than  thermo- 
dynamic precepts.  Here  it  is  applied  to  the  rapid  non- 
eguilibrium  expansion  of  atmospheric  water  vapor  in  air.  A 
method  is  established  for  demonstrating  the  validity  of  this 
theory,  originally  proposed  for  a mono-molecular  gas,  in 
the  case  of  a complicated  species  such  as  water  -.-ipor  , Cluster 
conditions  in  the  temperature  range  210-295®K  and  the  v.<»or 
pressure  are  calculated  throughout  the  collapse  of  a super- 
saturated metastable  vapor. 

An  effective  molecular  pair  interaction  energy  and  nearest 
neighbor  coordir.?i  ii>a  number  is  found  based  on  the  classical 
zero-poi  u enthalpy  of  sublimation  and  by  iteration  to  satisfy 
initial  equilibrium  li.mitioq  constraints  at  273°K.  The  internal 
energi'  redistribution  freqw.cy  results  from  close  matching 
of  the  experimental  pi  assure  distribution.  A method  of  deter- 
mining the  cluster  size  at  which  macroscopic  properties  become 
relevant  in  described;  the  size  is  found  to  be  a cluster  of 
about  122  molecules  at  273°K  based  on  a hard  sphere  model 

collision  cross  section.  For  the  water  vapor  dimer  an 

-21  -3 

-'dtilibrium  constant  of  4.0x10  (cm  ) is  found  at  273“K 
compared  to  3.1xlo”^^(cm~^)  obtained^^^  from  Keyes^^^  data. 

Two  models  for  water  vapor  clusters  result  (from  satis- 
fying all  constraints  imposed  on  the  theory)  which 
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are  given  in  terms  of  number  of  nearest  neighbors,  pair  inter- 


action energy  per  molecule  (2.9kcal/jj^^j^  at  273“iC)  and  possible 
structure  on  a cluster  by  cluster  basis  for  a classical  hard 
sphere  model  and  for  a model  approximating  a Pauling  type 
clathrate. 
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II.  FOREWARD 


The  work  described  here  is  an  attempt  at  developing  a 


method  of  applying  the  kinetic  theory  of  cluster  formation 
.(1) 


of  E.  R.  Buckle'  ' to  experiments  in  supersonic  nozzles.  The 
patient  support  of  the  Power  Branch  and  the  Atmospheric 
Sciences  Program,  Office  of  Naval  Research  over  many  years  and 
the  recent  support  of  the  Naval  Air  Systems  Command  is  acknow- 
ledged with  many  thanks.  This  effort  has  necessarily  been  an 
iterative  one  between  theoretical  development  and  numerical 
evaluation  over  more  than  ten  years  of  investigation.  Mathe- 
matical approximations  were  avoided  from  the  start  by  relying 
on  the  dio^tal  computer  from  the  start  of  development.  Time- 
shared  computer  facilities  were  provided  by  the  Naval  Weapons 
Lab,  Dahlgren,  Virginia,  as  early  as  1965  and  since  1967  by  the 
Naval  Academy.  The  importance  of  interactive  time-shared  faci- 
lities which  allow  the  investigator  to  treat  the  computer 
directly  as  a laboratory  simulator  is  not  to  be  underestimated, 
for  without  easy  access  tc  such  a system  this  work  could  not 
have  been  completed. 

The  cluster  theory was  developed  for  a binary  mixture 
of  monatomic  gases  and  it  seems  a bit  presumptuous  to  apply  it 
tc  as  complicated  a system  as  water  vapor  in  air.  However, 
the  best  data  and  most  clearly  understood  nozzle  behavior  was 
for  such  a mixture  and  at  least  a useful  first  approximation 
was  hoped  for.  Moreover,  the  phenomenon  of  condensation  is 
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universal  and  it  was  hoped  that  the  solution  could  be  worked 
out  in  terms  of  some  such  general  manner,  largely  independent 
of  the  exact  form  of  inter-molecular  potentials.  The  end 
result  is  far  better  than  expected  and  it  appears  that 
an  average  pair  interaction  energy  can  be  used  to  describe  as 
complicated  a non-equilibrium  interaction  as  that  bonding  a 
water  molecule  in  a cluster. 

The  objective  of  this  study  is  the  development  of  a 
technique  of  applying  the  cluster  theory to  the  condensation 
of  atmospheric  water  vapor  in  air  while  expanding  in  a super- 
sonic nozzle  (wind  tunnel) . Overall  behavior  of  such  flows 
7 8 9 

is  well  understood  ' ' now  and  accurate  measurement  of 
experimental  properties  provides  a means  of  comparison  with 
theoretical  results.  It  is  attempted  here  to  describe  the 
course  of  the  homogeneous  condensation  in  a nozzle  from 
molecular  rather  than  thermodynamic  parameters.  Parametric 
studies  of  the  influence  of  the  adopted  molecular  parameters 
are  given  demonstrating  the  sensitivity  of  their  influence 
on  thermodyneunic  properties. 

The  general  schema  is  to  establish  a valid  equilibrium 

initial  value  for  the  calculations,  then  by  calculation  follow 

the  course  of  the  expansion  in  a nozzle  through  the  region 
7 8 9 

of  condensation  ' ' . Two  models  of  water  clusters  evolve, 
one  a hard  sphere  collision  model  and  the  other  a clathrate 
type  model  more  akin  to  water  vapor.  After  establishing  what 
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appears  to  be  a valid  initial  equilibrium  condition,  the  non- 
equilibrium  expansion  follows  with  a close  matching  of  the 
experimental  pressure  to  that  calculated  from  the  theory. 

One  finds  that  the  non-equilibrium  cluster  behavior  is  dominated 
by  the  behavior  of  the  internal  enerfy  redistribution  frequency. 

The  calculations  described  here  were  perforaed  on  a H-635 
computer  using  the  Dartmouth  Time-Shared  System  (DTSS)  with 
the  long  non-equilibrium  calculations  being  performed  on  a 
Digital  Equipment  Corporation  PDP-15/40  medium  scale  computer. 
The  equilibrium  solution  requires  about  one  minute  while  the 
non-equilibrium  solution  requires  in  excess  of  one  hour. 

The  discussions  over  many  years  with  Dr.  W.  Dunning, 
University  of  Bristol  and  with  Dr.  P.  P.  Wegener,  Dr.  B.  T. 

Chu  and  other  Yale  University  Researchers  working  on  conden- 
sation, as  well  as  Dr.'s  B.  N.  Hale  and  P.  I>.  M.  Plummer  at 
the  University  of  Missouri-Rolla,  were  helpful  and  are 
acknowledged . 

This  work  could  not  have  been  completed  without  the 
continuing  advise  and  encouragement  of  E.  R.  Buckle,  Sheffield 
University  whose  help  was  sincerely  appreciated.. 

The  author  gives  thanks  to  the  Power  above  for  allowing 
a glimpse  into  his  very  nature  as  seen  in  the  dyneunic  beauty 
of  these  calculations. 
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III.  CLUSTER  DISTRIBUTIONS  AT  EQUILIBRIUM 


i 
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1.  Limiting  Constraints 

The  cluster  theory is  formulated  such  that  non-stationary 
growth  and  decay  rates  of  clusters  of  every  size  follow;  from 
phase  integral  specification  of  unimclecular  and  bimolecular 
reaction  probabilities.  Conceptually,  each  cluster  was 
originally  proposed  as  an  asscnbly  of  atoms  bound  together 
by  isotropic  London  forces.  Consequently  in  extending  this 
monomolecular  theory  to  the  polyatomic  water  molecule,  one 
must  accept  ein  average  value  of  the  pair  interaction  energy 
complementing  the  isotropic  London  forces.  As  will  be  seen 
later  this  does  yield  meemingful  results,  even  in  the  highly 
dynamic  non-equilibrium  situation  considered  later. 

Furthermore,  each  molecule  in  a cluster  can  act  as  an 
oscillator  with  the  number  of  vibrational  degrees  of  freedom 
per  cluster  Xg  determined  by 

Xg  = 3g  - 5,  for  the  dimer; 

Xg  = 3g  - 6,  for  all  larger  sizes; 

g,  being  the  number  of  molecules  in  a cluster.  The  vibrational 
energy  of  these  oscillators  circulates  within  a cluster  at 
frequency  v.  As  in  any  assembly  of  loosely  coupled  oscillators 
(with  nearly  similar  vibrational  characteristics) , at  any  lime 
a single  oscillator  may  possess  no  energy  or  alternatively 
all  of  the  energy.  Before  the  latter  condition  could  be 
realized,  however,  the  unimolecular  decay  criterion would 
cause  the  loss  of  that  particular  oscillator.  That  is,  when 
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the  energy  of  a single  surface  oscillator  exceeds  the 
"coupling"  energy  to  its  nearest  neighbors  it  is  lost  from 
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the  cluster;  this  energy  amounts  to  X u“,,  where  X is  the 

g 11  g 

nearest  neighbor  coordination  number  and  u“^,  is  the  pair 
interaction  energy. 

This  decay  energy  or  energy  increment  per  loss  (or  gain) 
should  not  be  confused  with  the  cluster  potential  energy  as 
in  addition  to  the  vibrational  energy  circulating  from 
oscillator  to  oscillator,  the  cluster  possesses  energy  binding 
it  together;  this  is  discussed  further  in  III-8  and  V-2. 

Another  parameter  to  be  specified  is  (0^)*  which  apportions 
some  fraction  of  the  Xg  overall  vibrational  modes  to  be  those 
of  surface  sites  eligible  for  decay.  The  limiting  constraints 
selected  to  obtain  this  parameter  as  well  as  those  above  are 
now  considered. 

The  cluster  theory  independent  parameters  to  be  specified 
or  determined  are; 

- pair  interaction  energy 

- cluster  internal  energy  circulation  frequency 

- surface  site  to  internal  site  parameter 

- nearest  neighbor  coordination  nianbers 

- monomer  concentration 

- cluster  size  of  effective  numerical  infinity 
*{  ) These  brackets  denote  the  set  over  all  g. 
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g>l 


a>l 


m\ 


where  C^  is  the  stationary  concentration  (f^  is  used  for  non- 
stationary states) . It  is  the  determination  of  these  quantities 
from  the  limiting  constraints  that  makes  the  solution  unique  to 
water.  (In  addition  the  collision  cross  sections,  molecular 
masses  involved,  and  some  knowledge  of  the  structure  is 
required,  see  II1-8) . This  is  in  marked  contrast  to  a nucle- 
ation  theory  developed  for  a specific  vapor (water)  where 
the  vibrational  characteristics  of  the  particular  molecule 
involved  are  included  as  parameters  uithin  the  theory. 

The  limiting  constraints  to  be  satisfied  here  for  saturated 
vapor  equilibrium  at  any  temperature  far  removed  from  the 


critical  point  are; 

(1)  Equilibrium  constant; 


» lim  o_ 


(2)  Zero  temperature  sublimation  enthalpy; 

(3)  Saturated  vapor  pressure;  p = Z C kT 

g>o  " 

(4)  Molecular  number  balance;  C,-  = Z gC 

g>o  ^ 

(3)  Incremental  surface  free  energy  change;  = 0 

g'+gto 

(6)  ‘Static  pressure  profile;  (pj^)  = {p..} 

theor.  exp. 

(calculated  by  starting  from  initial  stationary  (C^)) 

The  computational  procedure  described  next  is  iterative 
in  nature;  iterative  within  the  equilibrium  solution  as  well 
as  between  resulting  equilibrium  and  non-equilibritim  solutions. 
‘Non-equilibrium  condition  satisfied  later. 
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As  will  be  seen  later  in  V-2,  even  if  satisfactory 
molecular  descriptions  of  either  the  latent  heat  or  surface 
tension  were  known,  they  would  not  be  used  instead  of  limiting 
constraint  (6)  to  yield  a complete  evaluation  of  the  theory 
because  of  the  nature  of  v.  Approximate  means  of  calculating 
these  two  parameters  are  proposed  here  but  they  are  not  used 
as  limiting  constraints. 

The  first  of  the  variables  discussed  is  a^,  the  surface 
site  to  internal  site  parameter.  Recall  from  the  cluster 
theory that  this  parameter  apportions  the  fraction  of  total 
vibrational  degrees  of  freedom  that  are  decay  vibrations,  i.e., 
associated  with  surface  sites.  In  the  absence  of  a theoretical 
Interpretation  of  this  parameter  in  terms  of  the  eguations  of 
molecular  motion  within  the  cluster,  one  must  choose  a model 
for  the  cluster  structure.  The  very  nature  of  the  unimolecular 
decay  mechanism  in  Lite  theory  obscures  the  precise  effect  of 
the  internal  structure  of  clusters,  thus  detail  of  this  nature 
must  come  from  other  considerations. 


The  first  approach  is  to  envision  a water  vapor  cluster 
as  an  assembly  of  hard  spheres  with  collision  cross  sections 
defined  in  the  classical  sense^^^.  For  small  size  clusters 
all  molecules  would  appear  to  be  on  the  surface  (See  Figure 
VI-3.1)  and  paralleling  the  Pauling  type  clathrate  cage  where 
the  first  cage  closes  at  size  g = 20  we  define 

a = 1,  2%<20  III-l.l 


m 


The  cluster  theory yields  a limiting  value  of  for  clusters 
having  macroscopic  properties  in  the  equation  preceding  8.14^^\ 


a„  = lim  a = (-^)  III-1.2 

g->®  - ^ 

As  will  be  seen  later,  the  structure  based  on  hard  sphere  col- 
lision diameters  yields  a value  of  g^^  near  120.  When  = 120, 
= 0.98  so  the  parameter  remains  near  unity.  The  behavior 
imposed  here  for  what  might  be  called  the  hard  sphere  collision 
(HSC)  is; 


2 < g < 20 


Cg  = 1-0.1  (^),  g > 20 


III-1.3 


The  variation  of  a^,  g > 20  does  not  have  a great  overall  effect 
on  the  cluster  distribution.  Calculations  were  performed  with 
various  forms  of  the  expression  for  and  gave  only  minor 
differences  in  the  final  result,  the  same  cannot  be  said  for 
the  parameter  1^,  however  (see  below) . 

In  the  interest  of  closer  approximating  the  structure  of 
water  clusters  a second  structure  was  considered  (see  following 
discussion  on  collision  dieuneters  and  Section  V-2) . Implicit 
in  this  second  model  is  the  departure  from  Van  der  Waal  type 
forces  but  with  a structure  and  bonding  more  a)cin  to  water. 
Calculations  with  the  model  now  designated  as  the  water 
clathrate  model  (HCM)  show  a shift  of  to  about  g = 200 
so  becomes  about  0.83.  For  this  case  ve  have 
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III-1.4 
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Og  = 1 - 0.2  20  < g < 200 

It  will  be  seen  that  both  models  yield  good  results  in  pre- 
dicting the  dyneimic  behavior  of  non-stationary  states  (Section 
V-2)  . 

For  the  hard  sphere  collision  (HSC)  corresponding  to 
equation  III-1.3,  the  hard  sphere  collisir-r.  diameter  is^^^ 

a.  = r,  (i+g^-''^)  III-1.5 

A#  g * 

A closer  approximation  to  actual  dimensions  taken  from  a scale 
water  clathrate  model  - HCM  would  be 

Oj  g = r^{l+2(g)^/^),  III-1.6 

which  is  used  with  described  by  equation  III-1.4.  The  value 
•8 

of  r^  used  is  1.4.10  cm  for  all  calculations.  A summary  on 
all  aspects  of  both  models  is  given  in  Section  V-3. 

2.  Iteration  Scheme 

Having  specified  the  independent  parameters  of  the  theory 
and  the  limiting  constraints  which  must  be  satisfied,  it  is  now 
shown  how  these  conditions  are  met.  Referring  to  the  iteration 
scheme  of  Table  III-2.1,  one  initiates  the  problem  by  selecting 
values  of  temperature  T,  pair  interaction  energy  u'j^j^  and 
vibrational  energy  circulation  frequency  v.  A starting  value 
of  u®j^j^  was  obtained  by  subtracting  3kT  (for  rotation  and 
translation  modes  and  ignoring  expansion  work)  from  the  latent 
heat  per  noleaule  at  273®K  and  equating  the  result  to 
Selecting,  arbitrarily,  = 3.0  gives  a starting  value  of 
u*,,  = 2.98  k cal/mole. 
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TABLE  III-2.1 
TTERATI<ai  SCHEME 


choose 

T “K. 

choose 

(limiting  constraint  2) 

(A) 

choose 

V 

calc. 

C.  (perfect  gas  at  saturation  pressure 
^ for  p = CjkT) 

(B) 

choose 

Cl  (since  p = ^ Cg  kT) 

calc. 

calc. 

Fg  from  o^,  (a)  . (Limiting  case  g-^“» 

d)(V,  u°ii>) 

calc. 

Fig, 6)  from  (b) , (u„  = KLF(g,e)  » 

(go  to  A) 

t 

l^coinpare 

w(u'‘ii,lj 

F(0)  & P(g,e)  if  r(a)  s'  F(b)  within  e 
modify  v liimiting  constraint  1) 

choose 

g^  (Section  III-3) 

(C) 

choose 

mj^  (slope  of  Ig  vs  g) 

calc. 

(Cg) 

calc. 

P 

V 

(go  to  B) 

*— calc. 

=10-  “ '^v"’  'm  <• 

modify  (limiting  constraints  3 i 

(go  to  C) 

•—plot 

C,  if  A?  s<  o at  g = g^,  modify  mj^ 

(limiting  constraint  5) 

Use  {C  }as  initial  condition  for  non- 
9 

equilibrium  nozzle  calculation.  Follow 
expansion  down  to  onset  of  condensation 

at  pressure  p^^,  if  {p^)o,eor.  ^Pfe^exper.' 

modify  V {limiting  constraint  6) 

c JJi -.1  iti  r ^ iT-ru..A.7. 


Choice  of  the  vibrational  energy  circulation  frequency 

raises  questions  on  its  relation  to  that  of  vibration  spectra 

(4) 

of  water  observed  experimentally.  For  water  vapor  there 
are  frequencies  associated  with  shifts  in  mean  molecular 
position  as  well  as  with  individual  molecular  vibrational 
modes;  the  former  range  from  about  10^^  to  10^^  sec  ^ while 
the  latter  range  upwards  from  10^^  sec~^  to  well  above  10^^ 
sec”^. 

The  cluster  theory assumes  random  distribution  of  the 
vibrational  energy  Xg^T  circulating  at  frequency  v over  Xg 
vibrational  modes.  In  the  classical  approximation 

V « kT/h  II1-2.1 

(where  h is  Planks  constant)  such  that  the  circulation  fre- 
quency cannot  exceed  the  frequency  of  any  one  vibration.  It 
is  assumed  that  the  vibrational  energy  circulation  frequency 
here  is  more  closely  related  to  water  spectral  frequencies 
associated  with  intermolecular  vibrations  (shifts  in  molecular 
position) ; this  assumption  is  supported  by  later  results 
(Section  V-2) . At  the  initial  temperature  of  interest,  273'K, 
equality  in  equation  III-2.1  gives 

V = 5.7x10^^  sec"^. 

Next,  the  approximate  manomer  concentration  is  calcu- 
lated from  the  perfect  gas  law  (p  = CkT)  first  taking  C = Cj^ 
and  a linear  least  squares  fit  of  vapor  pressure  data  for  the 
liquid  state^^^.  It  is  assvuaed  then  that 
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p H £ C kT,  III-2.2 

g>o  " 

such  that  all  clusters  behave  as  a gas.  The  value  of  Cj^ 
calculated  above  fron  the  perfect  gas  law  must  be  reduced  to 
account  for  the  partial  pressure  of  the  clusters  and  at  the 
same  time  maintain  conservation  of  mass  by  molecular  nunjber 
balance  in  limiting  constraint  (4);  this  latter  condition 
requires  that  the  nvmiber  of  molecules  of  water  in  any  dis- 
tribution does  not  exceed  the  nurtoer  calculated  originally  from 
the  vapor  pressure  and  the  ideal  gas  law. 

Once  the  initial  value  of  is  chosen  at  the  start  of 

this  iteration  process,  the  value  of  1^  is  fixed  by  limiting 
constraint  (2).  The  next  step  follows  from  the  thfiery*^^  as 
the  equilibrium  constant  at  saturation  for  the  limiting  case 
(g->g<o)  is  = 1/Cj  . This  result  comes  from  the  definition 
of  the  parameter  J, 

J = Zn  Cj/u„,  III-2.3 

which  differs  slightly  from  In  of  the  saturation  ratio  because 
equation  II1-2.2.  At  saturation  J = 0 and  C,  = l/u_  = C,  so 
an  alternate  expression  for  J is 

J = Zn  Cj/Cjg.  III-2.4 

Perhaps  the  most  indirect  of  all  the  limiting  constraints 
specified  is  condition  (1)  but  this  gives  a meems  of  selecting 
a value  for  the  internal  energy  circulation  frequency  v that 
satisfies  other  conditions  of  the  theory.  Consider  first  the 
equilibrium  constant  which  results  from  detailed  balance  for 


each  g size  (see  {8},  reference  1 for  coinplete  details) 

u = C^Cj^Cg_j^  = KLgF(g,0),  III-2.5 

where  M,  P,  G are  convenient  grotipings  defined  by 


M=  (8im»j^j/m^)^/^rj^/v, 

111-2.6 

/x  -i\//x 

b(g>3,o,  = 

111-2.7 

\ fl  \ 

K ^ «'l,g/^l)^/(>'l.g/”'l)^'" 

111-2.8 

F(g,0)  = (e®/0^/^)  G(g,0), 

111-2.9 

G(g>3,6)  = {exp(e(Xg_j))-exp[e/(Xg-Xg_j-l)D 

/(1+x  e + x^  e^/,),  111-2.10 

where  0 = is  molecular  mass  of  water,  |j  is  the 

reduced  mass  and  X is  the  nearest  neighbor  coordination  number. 

In  the  thejrmodynaraic  limit  of  size^^^  (g-*-">) , 
ag=(36m/g)^'^^  and  Xg=X„  giving  the  asymptotic  values  independent 
of  size  g. 

L„  = (1/3611)^''^,  111-2.11 

F(g,0)  = F(0)  = lexp(X„0)-l/ 

0^-^^  (l+X^0+X^^0^/2) , 111-2 . 12 

= Mb^Fe  111-2.13 

Now,  making  use  of  the  thermodynamic  limit  of  size  (g+»)  at 
the  condition  of  a saturated  vapor  (J°0)  gives,  by  equation 
111-2.3,  111-2.11  and  111-2.13 

F(0)  = 111-2.14 

Is  ® 
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allowing  equality  here  with  equation  III-2.12.  Having  already 


defined  values  of  and  in  this  iteration  sequence  some 

value  of  V will  satisfy  both  III-2.12  and  III-2.14. 


At  this  point  some  is^licit  conditions  of  the  theory 


must  be  satisfied.  The  surface  free  energy  in  terms  of  the 


equilibrium  constants  is 


C = -kT  in(u./uj. 
i=2 


1II-2.15 


As  (I)  tends  to  u (g-«»)  and  the  clusters  take  on  bulk  properties 


(at  saturated  equilibrium),  by  equation  III-2.15  the  increment 


of  surface  free  energy  Ac  = c - c„_i  must  approach  zero;  this 
requirement  becomes  limiting  constraint  (5) . Also,  by  equation 
III-2.5  the  value  of  Ac  is  positive  when  G(6)/G(g,0) >Lg/L„. 

This  important  condition  begins  to  limit  the  range  of 
u®j^T  and  V.  In  general,  high  values  of  require  high 

values  of  v in  order  that  Ac  be  positive  over  all  g;  too  low 


a value  of  v causes  a hump  in  the  distribution  {C  } with 

g 

— 12  -1 

Ac<o.  At  the  upper  bound  on  v = 5.7*10  sec  , must 


be  reduced  to  about  2.91kcal/mole  to  satisfy  the  condition 


AC>o,  g-«»  as  well  as  equations  III-2.12  and  III-2.14. 


3.  Numerical  Infinity 


In  any  numerical  iteration  technique  an  upper  bound  on 


the  calculations  must  be  specified.  Selecting  the  value  of 


cluster  size  corresponding  to  numerical  infinity  is  considered 


next  although  tliis  in  reality  is  a separate  iteration  euid  an 


independent  method  of  determining  the  size  at  which  clusters 
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take  on  bulk  properties.  A first  approximation  at  the  value 
is  arrived  at  by  the  method  described  next,  further  refine- 
ment is  a result  of  iteration  to  satisfy  limiting  constraint  (5) . 

We  must  borrow  from  the  following  section  the  description 
used  for  surface  coordination  nusi>er  for  nearest  neighbors, 

Ig.  The  simplest  behavior  of  this  parameter  would  be  a linear 
variation  with  cluster  size.  The  slope  of  such  a variation 
is  taken  as  m^^,  which  is,  in  fact,  the  incremental  change  in 
surface  coordination  nuni>er  per  unit  increase  of  cluster  size 
(above  the  trimer) . The  equilibrium  computer  program  (See 
VIII -2)  is  now  run  with  all  variables  held  constant  except  for 
m^^.  An  arbitrary  value  of  this  parameter  yields  a minimum 
(artificial)  in  a plot  of  the  resulting  saturated  equilibrium 
cluster  concentration  versus  size  g.  Defining  the  cluster 
size  at  the  minimum  in  concentration  as  g=g*  and  plotting  this 
against  the  assumed  value  of  m^^  gives  the  result  in  Figure 
III-3.1.  The  difference  in  the  quantity  is  also 

plotted;  Xg,  being  the  value  of  Xg  corresponding  to  the 
minimum  in  {C  }. 

g 

As  the  value  of  slope  m^  is  decreased  the  cluster  size 
g*  at  the  minimum  in  {Cg}  tends  upward  in  size,  the  difference 
X^-Xg*  approaches  zero  and  as  is  shown  next  the  increment  in 
Ac  approaches  zero.  Below  a slope  of  m^  ° 0.01,  the  minimum 
in  g*  rapidly  exceeds  10^  (HSC) . Once  this  approximate  value 
of  is  determined  it  is  further  refined  by  iteration  to 
satisfy  limiting  constraint  (5). 
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In  Figure  III-3.1  the  effect  of  determining  g*  at  several 


-‘"V' 


I 


1 

I 


values  of  circulation  frequency  v and  is  evident  with  no 
great  change  occurring  with  small  variations  in  either  parameter. 
The  effect  of  slight  changes  of  slope  on  surface  free  energy 
are  shown  in  Figure  111-3.2  indicating  its  sensitivity  to 
structural  variation.  At  a cluster  size  of  about  g=123/  Ac 
is  effectively _ zero  (to  three  figure  accuracy)  for  the  hard 
sphere  collision  for  the  conditions  shown. 

4.  Nearest  Neighbor  Coordination  Kunber 

With  regard  to  the  variation  of  nearest  neighbor  coordina- 
tion number  only  the  behavior  with  linear  variation  is 
present  here;  numerous  other  descriptions  were  attempted  but 
satisfaction  of  the  limiting  constraints  was  not  achieved. 
Considering  the  previous  discussion  on  we  see  that  no 
matter  what  model  is  considered  (HSC  or  WCM)  the  first  20 
molecules  are  prescribed  as  surface  molecules;  the  immediate 
consequence  is  that  is  a surface  coordination  nueber  for 
g = 20.  Since  we  use  equation  111-1.3  to  describe  (HSC), 
g>20  nearly  all  of  the  other  sites  are  surface  sites  and  X^ 
must  be  essentially  a surface  coordination  number  throughout. 

It  will  be  seen  later  (V-2) , that  the  principal  contribution 
to  the  reaction  dynamics  is  from  small  cluster  sizes  (g<10) 
so  the  distinction  between  surface  euid  bul)c  coordination  is 
unimportant  for  either  the  HSC  or  WCM  case.  In  the  absence 
of  any  further  theoretical  relation  between  X„  and  a_  we 

g g 

assume  the  following.  The  upper  limit  of  the  set  {X^)  is  X^ 


js 

■V 

s 

S 
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whose  value  is  known  from  the  asymptotic  approximation  dis- 
cussed earlier  (equation  III-2.il  & I1I-2.13).  A lower  bound 
is  also  known  since  for  the  dimer,  aj  = 1 and  for  the  trimer 
Xj  = 2.  Thus  for  g>3  can  vary  only  between  the  limits  2 
and  X^.  For  water  (or  ice)  it  is  reasonab)e  to  suppose  X^ 

cannot  be  far  removed  from  a value  of  4.0  representative  of 

(4) 

tetrahedral  coordination 

Of  the  many  variations  of  X^  considered  the  following  two 
linear  models  were  most  successful  in  satisfying  the 
constraint: 

Xg  = Xj  + raj^(g-3);  mj^  = const.  ,3<g<g„, 

Xg  = X„;  g>g„,  1II-4.1 

and 

^g  = >-3  ^ 

Xg  = X^;  g>g„  III-4.2 

The  second  linear  case,  equation  III-4.2,  yields  results  quite 
similar  to  those  of  equation  II1-4.1  except  that  fine  adjust- 
ment in  m^^  allows  attaining  the  zero  change  in  incremental 
surface  free  energy  required  by  limiting  constraint  (5),  (See 
Figures  III-3.2,  and  III-4.1  to  3)  while  adjusting  g^  in 
equation  1II-4.2  does  not.  The  behavior  adopted  henceforth 
is  according  to  equation  III-4.1. 

The  final  limiting  constraint,  (6) , which  requires  agreement 
between  theoretical  and  experimental  supersonic  nozzle  static 
pressure  distribution  is  satisfied  in  section  V-2. 


LOG 


Figure  III-4.2 

Concentrations  at  Saturation,  273“  K 


Solution  of  the  equations  for  saturated  equilibrium 
(8.1  to  8.15,  ref.  1)  for  the  case  of  water  vapor  in  air  is 
given  in  a FORTRAN  interactive  computer  program  in  Appendix  2 
with  the  nomenclature  given  in  Appendix  1.  Corrections  to 
reference  1 are  given  in  Appendix  4.  A summary  of  HSC  results 
for  several  values  of  pair  interaction  energy  is  given  in 


Teible  I1I-5.1.  Results  are  also  given  for  the  dimer  equilibrium 
constant  for  comparison  with  that  obtained from  Keyes 


from  Keyes 


5.  Results  at  Saturated  Equilibrium 


When  one  succeeds  in  satisfying  limiting  constraints  (D 
through  (5),  the  results  appear  (HSC)  as  in  Figures  III-4.1 
to  III-4.3.  At  conditions  of  saturated  equilibrium,  satis- 
fying boundairy  condition  (5)  also  implies  that  AC^  = 0 as 
This  is  evident  from  the  equilibrium  distribution 

Cg  = Cj  expt(g-l)J-c^AT],  Ill-S.i 


which  at  saturation  (J=0)  is 


C = Cj  expl-C  ATI . 


III-5.2 


i.e.,  at  a given  temperature  the  variation  in  {C^}  at  saturation 
is  entirely  dependent  on  the  incremental  change  in  surface  free 
energy  AC^  with  cluster  size.  The  use  of  the  concept  of  an 
effective  numerical  infinity  is  tantamount  to  the  assur^tion 
that  clusters  of  size  g = g„  exhibit  bulk  properties. 

In  Figure  III-3.2  it  is  evident  that  the  distribution 
{Cg)  is  quite  sensitive  to  the  variation  in  (Ig).  A zero 
value  of  the  slope  roj^  in  equation  III-4.1  yields  a nearly 
linear  decreasing  distribution  with,  of  course,  no 

minimum  appearing.  A finite  value  of  m^  causes  a minimum  in 
the  distribution  with  limiting  constraint  (5)  requiring 

that  this  minimum  occur  at  effective  numerical  infinity.  In 
the  KCM  case.  Figure  III-5.1  shows  that  the  minimum  in  {C^} 
or  alternatively  Ac-*0  has  shifted  outward  to  the  neighborhood 
of  g„S200. 
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6.  Variations  with  Temperature 

The  computer  program  of  Appendix  2 gives  results  at 
saturated  equilibrium  at  any  temperature  by  varying  the  input 
temperature.  Having  satisfied  boundary  condition  (1)  through 
(5)  at  273°K,  the  parameters  listed  in  Table  III-S.l(c)  were 
used  at  the  other  temperatures  given.  The  monomer  concen- 
tration satisfying  boundary  conditions  (3)  and  (4)  was  the 
only  iteration  performed  with  the  end  result  that  condition 
(3)  was  satisfied  at  the  temperatures  given.  The  HSC  con- 
centrations obtained  are  given  in  Figure  III-6.1  with  the 
net  result  that  boundary  condition  (5)  was  satisfied  over  the 
range  290‘’K-230*k.  (See  later  discussion  on  v,  however,  in 
Section  V-2  with  regard  to  condition  (1)). 

7.  Surface  Tension  of  Clusters 

Theory gives  for  the  small  cluster  property  corresponding 
to  the  specific  tension  in  the  vapor-liquid  interface 

r = (?g/ga)  III-7.1 

which  differs  from  the  classical  result  obtained  by  equating 
the  product  of  surface  area  and  surface  tension  equal  to 
surface  free  energy,  i.e.j 

r = / (36r)^‘'^,  III-7.2 

where  in  effect  a = (36t/g;^'^^,  but  note  that  this  is  a for  a 
flat  surface  (cf.  line  3,  pg.  1280,  Ref.  1).  Eeurly  attempts 
at  evaluating  the  theory  included  satisfying  equation  III-7.1 
as  a boundary  condition.  Examination  of  results  in  Figure 
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CASE  (C 

.TABLE  III -5.1 

mm 

C,(CM'’) 

290 

4.80-  10 

273 

U60-  lO" 

250 

2.75 -lO'* 

230 

4.54  - 10  ” 

CLUSTER  SIZE  G,  HSC 


Figure  III-6.1 

Concentrations  at  Saturation,  290°  to  230°  K 


Fxio; 


CLUSTER  SIZE  G,  HSC 


Figure  III-7-1 

Surface  Tension  Variation,  Saturation  273°  K 


I 


I 
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III-7.1,  however,  shows  th^-t  for  the  Hard  Sphere  Collison, 
macroscopic  levels  7‘i  dynes/cm  at  273°K)  are  not  reached 
until  the  cluster  sire  approaches  g = 200  using  equation 
III-7.1,  or  much  larger  sizes  by  equation  III-7.2.  To  conform 
with  conditions  specified  thus  far,  for  g<g„  ^ 124  equation 
III-7.1  is  appropriate  and  for  larger  size,  equation  III-7.2 
would  seem  more  appropriate.  Results  in  Figure  III-7.2  for 
WCM  indicate  that  macroscopic  levels  of  surface  tension  would 
not  be  attained  in  either  case  until  sizes  in  excess  of  ctbout 
250  were  reached. 

The  net  result  is  that  the  use  of  stirface  tension  as  a 
limiting  oonstraint  precludes  the  possibility  of  satisfying 
all  other  limiting  conditions  imposed.  It  appears  that  the 
simple  molecular  viewpoint  used  to  interpret  the  equilibrium 
theory  here  requires  that  the  property  r corresponding  to  small 
clusters  be  considered  higher  than  the  macroscopic,  even  at 
the  cluster  size  satisfying  other  limiting  constraint  for 

8.  Latent  Heat 

As  a possible  alternative  limiting  constraint  the  thermo- 
dynamic property  known  as  latent  heat  was  considered.  In  the 
absence  of  a molecular  description  of  this  property  the  following 
crude  model  is  proposed. 

Consider  the  energy  added  to  a system  in  going  from  the 
liquid  to  vapor  state  (at  constant  pressure,  temperature)  as 
the  energy  used  in  breaking  bonds  in  the  liquid  to  establish 
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and  energize  an  equilibriint  population  of  nonomer,  dimer, 
trimer,  etc.  in  the  vapor  plus  the  work  in  moving  the  system 

boundaries.  In  the  liquid  state  some  of  the  nearest  neigh- 

(4) 

bors  may  be  hydrogen  bonded  in  tetrahydral  coordination  , 
others  may  not  be  involved  in  hydrogen  bonds.  The  fraction 

not  involved  in  such  bonding  nay  be  as  high  as  about  0.6  at 
(41 

273°K  with  most  results  predicting  a value  near  0.5.  A 

crude  attempt  is  made  below  to  allow  for  the  overall  effect 

of  such  unbonded  water  molecules  in  the  liquid  state. 

In  proposing  this  crude  model  of  latent  heat  some 

knowledge  of  the  structure  of  water  clusters  is  required. 

Cluster  structure  which  evolved  in  the  course  of  the  iterations 

to  satisfy  limiting  condition  (6)  is  discussed  in  section  V-4 

and  the  resulting  number  of  bonds  per  cluster  is  given 

in  Table  III-8.1  for  both  the  HSC  and  WCM.  The  bonding 

scheme  is  important  because  the  energy  stored  in  these  bonds 

must  be  considered  as  cluster  internal  potential  energy  in 

addition  to  its  vibrational  energy. 

TABI£  III-8.1 
BONDING  ARRANGEMCIT 


Liquid 

n ” ^00  = bonded 
n < 4,  unbonded 


! n,HSC 

n,wcM  1 

1 ° 

1 1 

? 

1 2 

i 

! 

1 

0 

1 

2 

g-1 

For  hydrcgen  bonded  molecules  in  the  liquid  state  with 
tetrahydral  coordination,  removal  of  a single  molecule  is 
taken  here  as  the  breaking  of  = 4 bonds;  in  the  non-hydrogen 
bonded  case  the  number  is  lower  and  is  discussed  below. 

Using  these  concepts,  for  each  dimer  formed  two  molecules 
have  their  bonding  changed  from  n=Xj^{liquid)  to  Oj  (vapor). 
Assuming  pairwise  additive  interactions  for  nearest  neighbors 
(at  absolute  zero  in  temperature)  the  energy  required  to 
separate  infinitely  a single  molecule  from  the  liquid  is 
X^u°j^j^,  In  forming  the  dimer,  two  such  actions  are  required 
and  the  energy  regained  on  bringing  two  infinitely  separated 
molecules  to  form  a dimer  is  Hence,  the  energy 

required  to  produce  a dimer  from  the  liquid  (at  T=0)  is 
(2X^-n2>u° By  similar  reasoning,  we  have  a g-mer: 

(aX  -n  )u°,,.  There  does  not  appear  to  be  any  reliable 
method  of  predicting  what  fraction  of  the  liquid  s*-ate 
molecules  are  hydrogen  bonded  in  tetrahydral  coordination, 
so  simple  reduction  of  X^( liquid)  from  4.0  to  a lower  value 
determined  by  iteration  is  introduced  to  give  an  overall  or 
effective  X^  in  the  liquid  state.  Thus  a resulting  value  of 
X^=  4.0  would  imply  100%  bonding  in  tetra-hydral  coordination; 
a lesser  value  implies  t)ie  presence  of  unbonded  molecules  in 
the  liquid  state.  Next,  to  adjust  for  temperature  above  zero 
we  must  add  an  amount  of  energy  equal  to  the  difference  in 
thermal  energies  of  the  g molecules  in  the  gaseous  cind  liquid 
states  as  well  as  the  work  of  moving  the  system  boundary 
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outward  from  the  liquid  state  to  that  encompassing  both 
liquid  and  vapor.  This  additional  volume  is  that  required 
to  allow  full  thermal  energy  for  each  cluster  in  the  vapor. 

We  can  now  express  the  T=0  energy  difference  from  liquid 
to  vapor,  say  for  the  dimer  as  (in  ergs  per  cc  of  dimer) : 

_ 1..0  bonds  bro)cen  (net)  ..  _ dimers  formed 
CgA.-njJu  11  aimer  fon^d * S Z 3 


III-8.1 


For  any  other  size,  including  the  monomer  we  have  (g^^'hg) 
CgU°j^j^  or  in  establishing  the  entire  equilibrium  population: 

III-8.2 


L»  = n,°  = £ (g>.„-n  )C  u«  , 

g>o  ^ g>o  y y 


= Xg  XT 


III-8.3 


where  Xg  ~ 3g-5  for  the  direr  and  3g-6  for  g>2.  At  constant 
pressure  and  temperature,  the  expansion  uock  pet  molecule  is 


where  Hg  = 5"^  WCM. 

In  addition  to  this,  energy  added  to  a system  to  vaporize 
it  must  do  vork  against  its  boundaries  and  increase  the  kinetic 
energy  from  that  in  the  liquid  state  to  that  of  the  vapor: 
here  the  kinetic  energy  of  the  liquid  state  is  neglected. 

Ac  in  the  theory the  clusters  have  vibrational,  rotational 
and  translational  energy  with  the  rotational  and  translational 
degrees  of  freedom  for  the  water  taken  as  6,  the  dimer  as  5 
and  all  others  are  taken  as  6.  The  vibrational  energy  is 
described  by^^^ 


h 

-'k 
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Thus: 

L'  = L°  + Energy {trans,  rotat.,vib.)  + Work 


L‘  = + I kT(C,)  -r  I kT{C2) 


6 . 


+ T £ C„  + kT  E x„C, 


g=2 


g=2 


3 3 


+ kT  E C_. 
g>o  ^ 


III-8.4 


To  comoare  this  result  with  accepted  data  in  the  literature. 


such  as  L = 535.5  cal/g  at  273°K^^^  , the  total  mass  ii„  in  grams 


for  the  population  {C^}  is  required.  Here 


3 


III-8.5 


where  is  the  mass  of  a single  molecule  and  V is  the  total 


volume  occupied.  Consequently  the  density  is  given  by 


V^/V  g/cm  . Thus  tabular  values 


(5) 


of  latent  heat  L are 


converted  to  cal/cm  by 


= fT, 


(L  cal/gram)  (ro^  1 3 C 9/^^)  • 


g>o 


III-8.6 


How,  equating  the  results  in  equations  III-8.4  and  III-8.6 


gxves 


u‘  , E (3^„-n„)C„  , '-1 

I,  = ^ g>o  S + I g 


5 kT 


g>o 


X g>o 


geo 

E c„ 


2 m.  E 3 C 
^ g>o  ^ 


^ 6 kT  q=3 


2 °1  ^ 3 Cg 

g>0  ^ 


!sl  q=2  ^ kT  q>o 

9 c m^  E g C 

g>o  ^ g>o  ^ 
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which  allows  direct  comparison  with  known  values  of  latent 
heat  in  units  of  cal/g. 

Iteration  for  X^(effective)  in  equation  II1-8.7  for  the 
liquid  state  at  273°K  yields  a value  of  3.04  with  L = 595.7 
cal/g  (H3C)  and  the  term  by  term  contribution  to  c is  given 
in  Table  III-6.2. 

TABLE  III-6.2 

Results  of  equation  111-8,7  ah  273*K,  L = 595.7  cal/g,  HSC  and  KMC 

Term  1 2 3 4 5 6 

Magnitude,  cal/g  478  88.2  0.004?  0.002  0.02  29.4 

% contribution  80.2  14.8  nil  nil  nil  4.9 

with  the  monomer  accounting  for  nearly  all  the  energy.  The 
solution  of  case  (c)  Tcible  111-5.1  was  employed  here  with 
identical  results  for  both  KCH  and  HSC  because  of  the  dominance 
of  the  monomer. 

It  does  not  appear  fruitful  to  speculate  further  on  the 
possible  fraction  of  unbonded  molecules  in  the  liquid  state. 
Reduction  of  1^  from  4.0  to  3.04  with  the  nearest  neighbor 
concept  employed  here  does  imply  unbonded  molecules  in  the 
liquid  but  the  fraction  is  still  unknown.  The  majority  still 
could  be  in  tetrahydral  coordination  with  an  apprecictble 
fraction  involved  in  fewer  bonds. 

Due  to  the  uncertainties  involved,  the  latent  heat  was 
not  used  as  a limiting  constraint  in  this  work. 
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IV.  CLUSTER  DISTRIBUTICWS  FOR  SUPERSATURATED 
AND  UNDERSATURATED  CONDITIONS 

1.  Experimental  Realizaticn 


The  system  considered  is  that  of  a dilute  solution  of 
water  vapor  in  a carrier  gas  of  air.  To  examine  the  states 
traversed  by  this  mixture  while  it  expands  adiabatically  in 
a nozzle  it  is  convenient  to  refer  to  Figure  IV-1.1  in  coordi- 
nates of  log  . and  log  T.  For  the  adiabatic  expansion  of  the 
mixture  p = cp^i  where  = Cp/C^  and  is  to  a close  approxi- 
mation equal  CO  the  ratio  of  specific  heats  of  air,  euid 
p = pRt.  The  slope  of  the  adiabat  (ABC)  is  only  a function 
of  the  specific  heat  ratio  (y^)  of  the  carrier  gas 

d (log  p)  3 T1 

a (leg  t)  IV-1.1 

while  fer  the  '/apor  the  saturation  line  (A*C')  is  obtained  from 
data^^^  or  the  Clausis-Clapeyron  relation  and  p = prt: 

log  Pj.  = 9.3286  - 2367.713Ar  300  < T < 250  IV-1.2 

The  vapor  and  carrier  gas  temperature  are  assumed  identical 
at  all  times For  a mixture  initially  at  point  A,  the  vapor 
state  is  undersaturated  and  J<0^^^  4.1 — 


Adiabatic  cooling  of  the 


mixture  occurs  on  expansion,  with  saturation  at  point  B. 
Beyond  point  (B)  the  vapor  is  supersaturated  witli  apparent 
collapse  of  supersaturation  beginning  at,  say,  C.  Cluster 
distributions  corresponding  to  thermodynamic  equilibrium 


from  point  A to  C are  now  considered  with  those  beyond  examined 


Point  A is  typical  of  the  initial  condition  encountered 
in  the  experiments  described  later.  The  undersaturated  con- 
dition at  point  (A)  where  J<0  follows  from  equation  IH-5.1. 
Since  J=ln  = In  the  monomer  concentration 

at  saturation  (A')  (at  the  temperature  of  point  (A)  must  be 
calculated  to  define  Conditions  at  (A*)  are  known  from 

the  solution  of  Appendix  (a)  with  iterations  satisfying  the 
saturation  pressure  equation  IV-1.2,  other  parameters  are 
retained  from  the  results  at  273®K.  Once  is  known,  the 
calculation  is  rerun  yielding  the  stable  undersaturated  (J<0) 
distribution  (a),  (296®K) , given  in  Figure  IV-1.2. 

For  the  supersaturated  condition  (J>0)  again  as  in  the 
undersaturated  case  the  equilbrium  equations  must  first  be 
solved  for  Once  is  known  the  solution  proceeds  as 

before  with  limiting  constraints  (3)  replaced  by  the  iteration 
for  the  isentropic  pressure 

= 0.452  (273)^1"^  IV-1.3 

Typical  results  for  temperatures  encountered  in  later  experi- 
ments are  given  in  Figure  IV-1.2.  As  noted  in  the  theory^^^ 
the  minimum  present  in  {C^}  when  J>0  corresponds  to  Volmer's 
critical  size  at  g=g*.  For  further  discussion  of  this  point 
see  Section  VII. 
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V.  NOS-EQUILIBRIUM  AHD  RELAXATION  CONSIDERATIONS 
. Constant  Volvune  Quenchin 


In  order  to  begin  to  compare  the  results  of  numerical 
analysis  with  experimental  data,  initial  equilibrium  conditions 
for  the  calculation  must  be  known.  Once  the  basic  predictions 
of  the  theory  of  equilibrium  etre  in  hand  it  is  possible  to 
eiq>loy  these  results  as  initial  conditions  and  begin  the 
calculations  involving  a gradual  departure  from  equilibrium; 
such  is  the  case  in  the  rapid  expansion  of  a condensable 
vapor  in  an  inert  carrier  gas.  However,  inasmuch  as  an 
adiabatic  expansion  is  difficult  to  interpret  on  the  molecular 
scale,  it  is  useful  to  examine  first  the  constant  volume 
quenching  process  and  then  to  proceed  to  the  case  of  a steady 
expansion. 

Consider  a mixture  of  saturated  vapor  with  a cluster 
population  (C^}  and  inert  gas  initially  in  thermal  equilibrium 
subjected  to  infinitely  rapid  cooling  to  a supersaturated 
state  and  maintained  thereafter  at  constant  volume.  The 
clusters  present  initially  at  equilibrium  are  thus  subjected 


to  quenching  with  the  rates  of  cluster  growth  exceeding  the 
rates  of  cluster  decay  tending,  in  effect,  towards  condensation 
with  an  exothermal  rate  of  heat  production. 

To  calculate  this  heat  production  consider  the  gas  phase 
reaction  by  collision  of  an  cluster  and  a Aj^  monomer  to 

give  an  A^  cluster 

Vi  ^ 


43 


in  the  interval  during  which  the  temperature  T is  constant 


and  the  total  volume  is  unchanged.  The  kinetic  theory  of 
cluster  formation specifies  the  energy  release  associated 
with  this  reaction  as  the  product  of  the  nearest 

neighbor  surface  cocrdination  number  and  the  Van  der  Waals 
pair  interaction  energy.  The  objective  here  is  to  specify 
the  net  exothermal  heat  of  reaction  associated  with  all 
reactions,  Eq.  V-1.1  for  a cluster  population  stib- 

jected  to  a rapid  quench  (or  heat  bath) . 

Several  approaches  to  this  problem  were  considered 
over  the  years  with  the  most  promising  result  given  here. 

To  describe  the  energy  of  reaction  of  V-1.1  an  energy 
balance  is  first  taken  on  the  bimolecular  mechanism 

\-l  *1  = \ 

and  then  followed  with  the  unimolecular  process 

\ = Vi  ^1 

within  a fluid  element  of  constant  volume,  temperature  and 
concentsyation  . 

Consider  first  the  bimolecular  reaction 

Vi 

proceeding  forward  at  rate  R^.  For  every  cluster  formed 

at  rate  R^  from  an  A , cluster,  1 monomer  is  'consumed'  and 
g g-1 

X u", , energy  of  reaction  is  involved  oer  A formed.  A 

g 11  - mer 

general  energy  balance  on  the  element  is  (where  t — dc/dt) 


^in  ^ *■  liberated 


^out  ^ 


3^ 


0 

stored 


with  bimolecular  input  equal  to 


^in,bi 

li]^ 


9-1  9 9-1  9 1 9 9 9 


V-1.2 


per  cluster.  The  reaction  rate  expressions  from  the  theory 
are  summarised  in  Table  V-1.1.  Energy  e®  is  the  zero  point 
cluster  energy  equal  to  the  sum  of  all  pair  interaction 
energies  within  a cluster  at  zero  temperature.  This  energy 
can  be  viewed  as  cluster  internal  or  potential  energy  involving 

all  molecules,  both  bulk  and  surface,  of  each  cluster.  Energy 
T 

E is  the  thermal  energy  corrections  for  temperatures  above 
absolute  zero  contributed  by  the  thermally-excited  degrees  of 
freedom  previously  defined  (Section  XII-8) . Thus,  for  the 
particles  involved  in  this  reaction  we  have: 

Particle  Energy 


A*  G”£  ® 

g 99  int.  g,rot.  g transl 

’'g  + I ?g  + I Sg  kT 

N-r  Vi  = ^-1  ^-^*7  «g-i  ^ ? ^g-1 

: E^  = Xj^  kT  + J kT  + J Sj^  kT  V-1.3 

D 

Finally,  we  have  c = ^s  the  energy  liberated  in  the 

bimolecular  reaction^^^ . It  is  important  to  focus  attention 
on  the  individual  particles  involved  as  it  is  the  (g-1)  size 
cluster  and  the  monomer  which  enter  the  fluid  element  while 
the  reaction  "within"  the  element  is  written  for  the  (g)  size 
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cluster.  In  essence  we  make  an  energy  balance  (biroolecular 


part  of  reaction)  only  on  the  new  particles  df^  = R^dt 


entering  the  element  as  f^  remains  constant  with  the  element. 


The  energy  out  of  the  fluid  element  associated  with  the 
bimolecular  reaction  is 


£„...  = e°  r"^  - dq*^/dt, 

out,br  g g g g • 


V-1.4 


where  a positive  value  of  heat  input  is  in  accord  with  the 
first  law  of  thermodynamics  for  the  total  energy  of  the  gas: 


da 

3t 


al- 


dw 

at 


V-1.5 


Next  we  consider  the  unimolecular  reaction 


VS-1  "1 


within  the  element  proceeding  at  rate  and  involving  the  decay 


of  df  = R dt  g size  clusters.  The  initial  energy  before 
g»u  g 


reaction  is 


. 0 + T + 

c.  =e  R + € R„ 
in,u  g g g g 


V-1.6 


with  reaction  energy  e = 1 u°,,  we  have 

g g 11 
■g 


= -"a  * 

9 9 


So 


c.  = e°  R + c R^  + R , 
in,u  g g g g g g 
lib. , 


V-1.7 


The  energy — after  the  <jniroolecular  event — leaving  the  fluid 
element  is 


-I 


'final, u ~ ^g-1  *^9  * ^9”1  ^9 


V-1.8 


+ e£  Rg  - dq  /dt 

*having  taken  heat  liberated  within  the  element  as  positive 
for  bimolecular  reaction 
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Rewriting  the  general  energy  balance  on  the  fluid  element 


^®out,bi  ” *'in,bi^  ^ ^®final,u  ^initial, - 0 V 1.9 
and  subiStituting  Eqs.  V-1.2,  V-1.4,  V-1.7  cmd  V-1.8  gives 


+ (e°  Ra  ■ 

g g g g 

- (e°  , R^  + - R^  + e?  R^  + e®  R*) 

g-1  g g-i  g i g g g 

^ ^ ^ "g-i  \ ^ ">g  ■ 

- (e°  R^  + R^  - Rj 

g g g g g g 


= 0 


bi. 


Uni. 


So 


- dq‘’/dt  - e°  (R*  - r")  - e^_,  (R^  - r")  - eT  (R*  - R") 

g-1  g g g-1  g g 1 g g 


dq“/dt  - C«  (R+  - R-) 


g g 
= c 


Defining  dq  = dq^  + dq'^  gives 


H '^g  - - ^g-1  ««g 

- '^g  <<  - 


RJ  - e^  (R^  - r") 

g 1 g g 


v-i.io 


^ ~ *n 

Further,  if  (R  - R ) = A , or  net  flux  on  the  element,  Eq. 

g g g 


V-1.11 


V-1.10  can  be  written  as 

i = - <^^-1  - ^g-l  ^ «g 

where  e®  , = \ lU”,,  (HSC)  or  (g-2)u®„  (WCM) . Conceptually, 

y“A  £ S** 

the  balance  on  the  fluid  element  now  appears  as  in  Figure  V-1.1. 
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This  non-stationary  result  for  dq/dt  is  not  related  to 
the  latent  heat,  ctn  equilibrivim  constant  pressure,  temperature 


concept.  It  does  define  the  reaction  energy  on  net  growth 
or  decay  of  clusters  which  behave  as  a perfect  gas  and  is 
the  heat  addition  term  in  the  first  law,  Eq.  V-1.5.  In  a 
binary  system  such  as  water  vapor  in  air,  care  must  be  taken 
in  precisely  defining  the  overall  system  described  by  equation 
V-1.5  since  the  heat  production  described  by  equation  V-1.5 
is  considered  as  originating  in  an  external  source. 

If  is  the  specific  heat  per  cc  of  the  carrier  gas, 
the  specific  heat  of  the  mixture  is 


*=v  = ^vo  ^ ^ (S^  + Cg)/2) 

garj. 


V-1.12 


where  f denotes  a non-stationary  concentration. 

Combining  Eqs.  V-15  and  V-1.11  for  the  constant  volume 
cases  gives  for  the  mixture: 


V-1.13 


g "g 

If  p°  and  T®  are  the  pressure  and  temperature  before  quenching 
the  perfect  gas  law  is 


^ = El  ^ V-1  M 

9a> 

and  p = kT  (f^  + Z f^)  V-1.15 

Equations  V-1.13  and  V-1.14  are  non-in tegrable  so  that  p 
and  T must  be  solved  for  numerically  from  an  initial  equilibrixim 
cluster  distribution  {C^)  at  p°  and  T°.  The  pressure  may  be 
arrived  at  either  by  V-l.H  or  V-1.15  where 


V-1-1) . 


so 


(see  also  Table 


dt^=  % - ^ g+1  ^ « g+1  - g- 
With  this  set  of  equations  one  may  evaluate  in  a step-wise 
manner  the  changes  in  state  that  occur  in  a constant  volume 
containing  moist  air  as  it  is  quenched  or  heated.  The  tem- 
perature change  produced  by  the  entire  population  {f^}  shift 
is  applied  after  calculating  equation  V-1.16  for  all  g with 
the  same  instantaneous  initial  conditions. 

The  non-equilibrium  behavior  of  the  initial  stationary 
cluster  distribution  {C^}  described  by  case  c.  Table  III-5.1 
eund  subjected  to  an  instantaneous  temperature  drop  of  0.0C01*K 
(the  approximate  average  temperature  drop  per  collision 
interval  in  later  experiments)  is  given  in  Figure  V-1,2. 
Details  of  performing  this  calculation  are  given  in  the 
following  section  and  in  Appendix  (3);  all  expansion  terms 
are  simply  nulled  out  for  this  constant  volume  case  and  the 
output  is  examined  each  calculating  interval  of  0.<30  n sec. 

In  order  to  gain  some  insight  into  the  system  behavior  in 
the  simplest  case,  the  temperature  changes  indicated  by 
equation  V-1.13  (practically  nil)  are  also  suppressed  after 
the  initial  quench  of  O.OOOl'K. 

The  ordinite  of  Figure  V-1,2  is  obtained  by  dividing  each 
rate  df^  by  the  final  positive  value  calculated  for  that  size 
g (defined  df^)  while  the  abscissa  is  the  nuidSer  n of  tiine 
increments  At  = 0.4 80  n sec.  The  dimer  behavior  is  not  shown 


because  of  relative  magnitude;  at  n=l,  df-  = -1.22xiC';  at 


o 

n=40,  dfj  = -0.708x10  . The  diiiier,  tlierefore,  behaves  here 
much  like  a monomer  since  its  rate  is  negative  throughout; 
larger  size  clusters  all  have  rates  which  are  first  positive, 
then  negative.  It  will  be  seen  later  that  with  expansion  the 
dimer  rate  is  positive  until  departure  from  the  adiabat  occurs. 

On  chainging  AT  to  a positive  increment  of  +0.0001“K,  or 
a heat  bath,  the  opposite  behavior  is  observed  (though  not 
symmetric)  with  the  dimer  increasing  throughout.  The  dimer 
literally  receives  a "shower"  of  fractured  trimer  leading  to 
positive  rates.  If  AT  is  increased  slightly  to  say  0.001  or 
0.01,  then  a negative  rate  with  the  typical  type  of  behavior 
for  the  larger  cluster  sizes  shown  in  Figure  V-1.2.  The  time 
to  the  first  reversal  of  sign,  i.e..  at  is  given  for 

each  cluster  size  in  Figvire  V-1.3  as  a matter  of  reference. 

The  behavior  of  the  cluster  is  evident  in  Figure  V-1.2 
with  a wave-like  disturbance  oscillating  in  the  medium  until 
dcunped  by  the  growth  and  decay  processes.  Interpreting  the 
results  of  Figure  V-1.3  as  relaxation  times  (in  attaining 
equilibrixnn)  would  be  incorrect  since  gradients  would  still 
exist  in  the  medium  at  the  times  shown. 

A reassuring  check  on  the  theory  cind  computer  analysis 
is  available  at  this  point.  Since  the  equilibrium  solution 
was  obtained  by  detailed  balance  involving  non-equilibrium 
rates,  the  use  of  cin  equilibrium  distribution  {C^}  at  tem- 
perature T in  a non-equilibrium  calculation  for  {f^}  at 
temperature  T should  yield  stationary  rates  with  = r”^. 
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1 2 3 4 5 6 7 8 9 X)  11  12  13  14  15 

CLUSTER  SIZE  G.  HSC 

Figiite  V-1.5  „ 

Reversal  Tines  for  Quench  of  AT  ==  -0.0001  K 
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all  g.  The  use  of  hard  sphere  collision  diameters  gives  a 
check  to  three  digit  accuracy.  It  was  found,  however,  that 
a check  to  greater  accuracy  was  possible  by  modifying  the 
monomer,  dimer  collision  cross  sections,  calculating  the 
equilibrium  distribution  (C^},  then  checking  for  stationarity. 
For  the  hard  sphere  clathrate,  the  use  of  j = 1.95  r^ 

(vice  2.0  r^)  and  = 2.05  (vice  2.0  r^^)  gave  equality 

in  = R~g  from  three  to  five  digits,  all  g,  with  the  results 
at  273°K  given  in  Tcible  V-1.2. 

For  the  constant  volume  quench  (or  heat  bath)  it  would  be 
possible  to  follow  the  behavior  of  individual  clusters  in 
attaining  a new  equilibrium  by  continuing  the  calculation  of 
Figure  V-1.2.  If,  however,  a new  ten^serature  Icwer  (or  higher) 
by  another  increment  of  AT  is  imposed  at  time  t At,  it 
becomes  impossible  to  interpret  what  is  happening  on  an 
individual  cluster  versus  time  basis.  This  is  the  process 
occurring  in  a rapid  expansion  where  the  rate  processes  of 
each  time  increment  modify  those  of  the  succeeding  one.  It 
appears  fruitless  to  try  to  relate  such  results  to  relaxation 
times,  rather,  investigation  of  the  behavior  of  the  cluster 
population  on  the  whole  is  more  rewarding. 
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TABLE  V-1.2 

STATIONARY  GROWTH  AND  DECAY  RATES  AT  273‘>K 


5 R^(cm~^  sec~^) 

R {cm~®sec~^) 

2 0.40966  • 10^^ 

0.40958  • 10^^ 

3 0.26336  • 10^^ 

0.26337  • 10^^ 

4 0.47028  • 1C^° 

0.47020  • 10^° 

5 0.12430  • 10^® 

0.12429  • 10^® 

6 0.25205  • 10^® 

0.25201  • 10^® 

7 0.43719  • 10^^ 

0.43711  • 10^^ 

■ 8 0.68725  • 10® 

0.68712  • 10® 

9 0.10143  • lo’ 

0.10142  . lo’ 

1 10  0.14395  • 10^ 

0.14392  • 10^ 

Adiabatic  Expauision;  Comparison  with  Experiment 

The  remaining  limiting  constraint.  No.  6,  is  now  considered. 
In  the  expansion  of  moist  air  in  a supersonic  nozzle  pressure 
"humps’  which  are  steady  in  time  develop  in  the  regions  of 
condensation.  The  pressure  distribution  selected  here  for 

comparison  is  free  of  shock  wave  pressure  disturbances  which 

7 3 9 

could  be  within  a zone  of  condensation  ' ' . Consequently 
the  entire  difference  in  static  pressure  from  that  encountered 
in  a moisture-free  nozzle  flow  is  due  to  the  effects  of  heat 
addition  from  condensation.  He  know  therefore  that  a matching 


'll 
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of  calculated  pressure  with  the  experimentally  measured 
pressure  should  be  a meaningful  test  of  the  parameters  involved 
in  the  theory.  One  further  test  (which  is  the  subject  of 
future  investigation)  is  the  prediction  of  other  experimental 
results  based  on  the  parameters  dete2nnined  here. 

Having  gained  some  insight  into  possible  modes  of  cluster 
behavior  it  is  worthwhile  referring  to  Figure  rv-1.1  to  define 
the  approach  used  in  comparing  numerical  calculations  with 
experiment.  For  a given  set  of  molecular  parameters  it  is 
possible  now  to  calculate  conditions  at  saturation, 

9 

point  A',  then  calculate  the  undersaturated  initial  cluster 
distribution  both  at  initial  ten?)erature  ^01-  Using 

{Cg)j^g  as  an  initial  distribution  it  will  now  be  possible  to 
calculate  the  change  in  cluster  concentration  df^  for  all  g, 
by  rate  expression  V-1.16  modified  to  include  the  effects  of 
expansion.  The  partial  pressure  of  the  entire  population 


{fg}  rs  then  simply 


Pv  = gto  ^g 


The  objective  here,  therefore,  is  to  calculate  this  pressure, 
add  to  it  the  local  pressure  of  the  carrier  gas  and  compare 
the  result  with  the  static  pressure  obtained  from  experiment 
along  the  equivalent  path  .A-B-C,  Figure  IV-1.1.  To  do  this 
we  must  add  to  the  results  of  the  previous  section  the  effect 
of  expansion  occurring  in  the  flow. 

For  the  steady  inviscid  expansion  experiencing  no  body 


forces  the  net  rate  of  internal  energy  change  per  cc  is  given 
by  V-1.5  where 


V-2.2 


Considering  one  diriensional  steady  flow  at  velocity  u in  a 
nozzle  of  cross  sectional  area  A,  density  p,  the  mass  flow 
rate  obeys 


^ (pAu)  = 0,  V-2,3 

and  the  volumetric  flow  rate  is  simply  Au.  The  quantity 
d(ln  v)/dt  is  thus  given  by 

(In  v)  = (In  Au).  V-2.4 

Since  d In  y/dt  = 1/y  dy/dt,  hence 

In  Au  - ^ [A^  + u||j  = ^ [In  u + In  A]  V-2.5 


Now,  rewriting  equations  V-1.5,  V-2.2  and  V-2.4  and  defining 
de  = c^  dT  gives 


dr  _ ^ _ d(ln  Au) 
dt  ~ dt  " 3t  ' 


V-2.6 


and  having  dq/dt  by  equation  V-1.6,  c^  by  equation  V-1.12  and 
p by  equation  V-1.15  allows  expressing  the  temperature  in 
equation  V-2.6  by 


dT  _ 1 , - .0  . ,T  4.  ,-T  . _R,  in 

^ % ''g>l  9-1  ^ S-l  ^ "l  ^ ^ 


- kX  (f  + Z f ) 
° g>0  9 


d(ln  Au) 

— at — 


•1. 


V-2.7 


The  pressure  follows  from  the  ideal  gas  law  since 
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dp ddn  T/v, d 

at  - P 3t Pdt  ' 

so  by  equation  V-2.4 


V-2 . 8 


dp ,d(ln  T)  ^ Cn  Au) , 

dt  - Pt-dt ^ “dt 


V-2. 9 


or  combining  equations  V-2. 9 and  V-2.1  allows 


d In  p _ 1 , 

dt  “ C T 


g>l 


(e 


.1.  J.  J.  IIJ 

+ c , + e,  + e)6 

0-1  1 o' 


g-l  ^ '-g-1 


- [1  + ^ (f^  t E f ) 
V g>o 


V-2. 10 


Hera  again  equation  V-2. 7 and  V-2. 10  are  non-integrcdsle 
and  numerical  solution  is  necessary.  As  in  the  constant  volume 
case  the  pressure  may  be  found  from  equation  V-1.15  rather  than 
from  V-1.14. 

One  further  consideration  is  needed  to  complete  the  des- 
cription of  the  expanding  medium.  Since  the  volume  of  the 
system  varies  with  time,  a term 

V-2. 11 


- f (d/dt)  In  V 
9 


must  be  added  to  eq.  V-1.16,  while  for  the  diluent  gas  we  have 

V-2. 12 


- f^  (d/dt)  In  V, 


noting  that  the  volume  change  is  already  specified  in  terms 
of  the  area  and  local  flow  velocity  in  eq.  V-2. 5.  The  change 
in  concentration  of  any  cluster  size  during  expansion  is  thus 
written  as  the  sum  of  Eqs.  V-1.16  and  V-2. 11 


dt 


- f. 


(In  V)  . 


V-2. 13 


1 
i 

? I 


For  the  monomer  rate  equation,  we  have  (constant  volume) 


df,  = - E g(ffg)  at. 


at^’ 


V-2.14 


or  with  expansion 


9=2 

,df„,  ..  . a 


df.  = - Z q(“*-g)  at  - f,  ^ (In  V), 
^ g=2  at  ^ 


V-2.15 


Adiabatic  expansion  of  moist  air  in  a supersonic  nozzle 
can  now  be  examined.  The  case  considered  here  (Exp.  405, 
Ref.  8)  is  free  of  pressure  discontinuities  from  shock  waves 
so  the  only  deviation  of  static  pressure  from  that  of  a dry 
adiabatic  expansion  is  that  caused  by  the  heat  input  of  con- 
densing vapor.  The  information  required  of  the  experimental 
conditions  for  input  to  the  numerical  calculations  is  the 
initial  pressure,  temperature,  the  nozzle  profile,  i.e.,  the 


area  A^^^j  and  slope  dA/dx.  In  order  to  )aiow  the  time  per  data 


interval  the  flow  velocity  u^^j  in  the  nozzle  is  required; 

,(7,8) 


this  follows  from  solution  of  the  equations  of  motion 


(1) 


Inputs  from  the  equilibrium  theory  are  the  initial  cluster 


distribution  {C^}  and  the  accompanying  parameters  (Say  case  c. 


Table  III-5.1).  The  objective  is  to  match  the  pressure 
deviation  in  the  nozzle  observed  at  onset  of  condensation — 
thus  satisfying  limiting  constraint  (6) ; 

The  FORTRAN  computer  program  for  calculating  the  non- 
equilibrium rates  of  change  of  cluster  concentration,  the 
temperature,  pressure  and  other  parameters  is  described  and 
given  in  Appendix  3;  the  general  flow  ediart  is  also  given. 


i i 


sc 
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In  essence  the  program  calculates  the  rates  of  growth  and 
decay  for  each  cluster  size  from  the  dimer  upwards  in  each 
calculating  interval  thus  giving  the  concentration  of  each 
species  with  respect  to  rime  (or  location) . The  static  pressure 
reported  is  thus  the  sum  of  the  partial  pressures  of  all  species 
present,  from  the  monomer  upwards  and  the  carrier  gas. 

The  calculating  interval  used  requires  further  consideration. 
The  earlier  discussion  on  cluster  behavior  under  constemt  volume- 
quench  conditions  showed  that  time  increments  on  the  order  of 
tlie  collision  interval  were  required  in  that  case.  With  the 
complication  of  expansion  the  interval  to  be  used  was  uncertain; 
the  minimum  was  still  the  collision  interval  but  perhaps  a 
longer  one  could  be  used  to  conserve  computer  time. 

Collision  time  corresponding  to  the  mean  free  path  is 
calculable  from  the  collision  numbers  of  the  theory^^\  i.e.: 

^g  = '^'og  * ^'ig'”^'  tX^)"!  V-2.11 

where 


and 


f 

f 


0 

1 


Z 


x,g 


.8rlcT,l/2  2 


X = 0,1 


Using  concentrations  from  typical  calculations  we  find,  for 
example,  that  at  270“K,  “v  4nsec.  This  dictates  a scheme 
of  variable  step  size  which  is  described  further  in  Appendix  3. 


There  appears  to  be  no  significant  difference  in  results, 
however,  between  a constant  step  size  and  a variable  step 
size.  It  is  possible  by  this  method  of  variable  step  size, 
though,  to  reduce  computer  time  from  about  16  hours  per 
solution  to  about  1.5  hours.  A step  size  twice  the  largest 
collision  interval  will  lead  to  an  unstable  numerical  solution. 

A fxirther  area  requiring  clarification  is  the  experi- 
mental input  data.  For  the  case  considered  here  data  on  Area, 
Velocity,  etc.  were  for  0.1  cm  intervals.  This  j too  coarse 
and  must  be  refined;  simple  linear  relations  were  used  between 
data  intervals. 

A final  point  of  clarification  concerns  the  initial  cluster 
distribution.  Referring  back  to  Figure  IV-1.1,  moist  atmos- 
pheric air  entered  the  experimental  apparatus  in  condition  A. 
Some  expansion  already  occurred  at  the  point  where  experimental 
pressure  was  first  measured.  The  temperature  had  already 
dropped  to  273°K  at  this  point  from  295. 8°k.  To  account  for 
the  change  in  cluster  concentration  from  point  A to  B,  two 
approaches  were  tried. 

First,  the  undersaturated  distribution  at  A was  calculated 
and  expanded  adiabatically  to  273°K  considering  the  absolute 
number  of  each  species  to  remain  unchcinged  (frozen  equilibrium) . 
Second,  the  equilibrium  distribution  at  273°K  was  calculated. 
Inasmuch  as  the  first  method  gave  vapor  pressures  slightly 
lower  than  those  on  the  adiabat.  Figure  IV-l.l,  Uie  second 
condition  was  used  as  the  initial  condition  in  the  results 
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reported  here.  There  was  very  little  difference  in  the  two 
distributions  and  both  would  give  nearly  similar  final  results 
except  that  the  first  case  would  follow  a path  just  under  the 
adiabat. 

A bit  of  historical  development  is  injected  at  this  point 
to  show  how  the  final  solution  was  attained.  In  evaluating 
numerical  results,  initially,  the  greatest  confidence  was 
placed  in  the  results  of  tne  rate  equation,  V-1.12.  The 


argument  stemming  from  the  fact  that  and  closely 


related  through  detailed  balance  at  equilibrium;  with  most  of 


y.e  experience  gained  in  the  many  variations  of  fCg).  From 


this  the  mass  fraction  of  "condensate”  is  defined  by  considering 
the  evolution  of  all  clusters  from  the  dimer  upwards.  The  mass 


fraction  of  water  clusters  per  kg  of  air  is  cbtained  by 

-3 


•10 


9 


g>l  f_ 


V-2 . 12 


o ' o 

where  m^^  and  m^  are  the  mass  per  nonomer  and  mass  per  air 
molecule  (avg)  respectively;  the  initial  value  of  X = 0.0C6 
being  neglected.  If  this  exprsssion  is  suroned  to  include 
the  monomer,  the  mixing  ratio  'sometiro.»s  called  the  specific 
humidity)  is  obtained. 


g>o 


V-2. 13 

"b  ^ ” *o 

should  remain  constant  throughout  equal  to  the  initial  value, 
here  equal  to  4.86  g/kg.  This  is  in  effect  another  internal 
boundary  condition  to  be  satisfied. 


I i 
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The  basis  of  consparison  for  X is  the  result  of  solution 

7 8 

of  the  equations  of  piotion  for  the  flow  and  given  elsewhere  ' . 
Although  this  earlier  solution  relies  on  the  use  of  macroscopic 
latent  heat  and  therefore  has  limited  validity  with  respect 
to  the  condensate  mass  fraction,  it  does  appear  to  be  a useful 
first  approximation  and  is  plotted  for  experiment  #405  in 
Figure  V-2.1.  Some  early  results  for  X are  given  iii  Figure 


An  interesting  problem  is  encountered  at  this  point, 
namely,  the  strong  dependence  of  X on  u.  If  case  c,  HSC, 

Table  III-5.1  is  used,  curve  (c)  results  which  is  quite  dis- 
appointing in  that  little  or  no  condensate  appears;  the  small 
amount  appearing  is  late  in  tine  or  downstream  of  experimentally 
observing  results.  If  a temperature  dependence  is  imposed  on 


log  'J  = ^ log  5.55  ' 10^^, 


with  T^  = 273,  still  satisfying  case  (c)  Table  III-5.1,  then 
curve  IB)  results  which  is  premature  or  upstream  of  observed 


/ 


...  this  point  an  int.ense  investigation  was  made  on  the 
influence  of  other  parameters  in  the  theory  as  well  as  different 
initial  condition  (Table  II  1-5.1).  Referring  to  the  rate 
expression  in  Table  V-i.l  one  finds  that  for  curve  (c)  the 
decay  rates  are  too  large  and  the  system  evolution  is  incomplete 
by  the  ti.ee  the  nozzle  exit  at  1 = 9 cm  is  reached.  Farametric 
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Figure  V-2.1 

Condensate  Mass  Froction,  HSC 


studies  showed  there  was  little  possibility  that  variation 
of  Yt  ®tc.  could  substantially  affect  the  system 

evolution  so  attention  was  turned  to  v. 

It  appeared  now  that  the  frequency  v was  the  dominant 
parameter  involved.  Examination  of  Figure  V-2.2  shows  some 
of  the  variations  of  v that  were  considered,  with  curve  (c) 
and  (b)  corresponding  to  those  in  Figure  V-2,1.  Since  the 
experimental  results  appeared  to  be  bracketed  by  cvirves  (C) 
and  (B) , intermediate  variations  for  v were  considered  with 
the  following  yielding  the  most  promising  results  (curve  D) 


log  V = log  V - 


‘o  "R'^IO' 


V-2.15 


Having  gained  some  insight  into  the  nature  of  the  behavior 
of  the  theory's  pareuneters  attention  was  shifted  to  comparison 
of  numerical  results  with  the  actual  experimental  pressure 
distribution.  Figure  V-2.3  shows  the  degree  of  con^arison 
in  static  pressure  (dimensionless  P/Pg^)  when  n = 5/4  and 
Tj^  = 213‘’K,  keeping  in  mind  that  the  calculations  were  initiated 
at  T = 273®K  at  X = 1.5  cm.  The  isentrope  is  followed  numeric- 
ally as  is  the  experimentally  (^served  pressure  deviation. 

The  result  in  Figura  V-2.3  for  X was  expected  to  differ  from 
the  earlier  result  as  pointed  out  previously  because  of  the 
earlier  use  of  latent  heat. 


Once  the  close  match  in  pressure  was  achieved  reexamination 
of  parameters  in  the  theory  showed  nearly  direct  proportionality 


between  v and  collision  nxjrber  Zj^j^  for  the  monomer-nsonciaer 


interaction.  If  fact  allowing 

V = k Z 


11 


V-2 . 16 


gave  a match  in  static  pressures  within  a few  percent.  A 
small  variation  in  k (from  2.05x10“^^  to  1.21x10”^^)  with 
temperature  and  saturation  according  to 


k = - 7.431  - 0.0789  s + 9.54 

T 


V-2 . 17 


yielded  a pressure  distribution  nesurly  identical  with  experiment. 
(See  Figure  V2.3). 
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15  hours  of  actual  computer  time  (CPU)  on  a medium  scale 
computer  (POP  15/40),  This  is  reduced  to  about  1.5  hours  by 
using  variable  collision  interval  basis  described  eeirlier 
and  in  more  detail  in  Appendix  3.  As  any  computed  parameter 
is  observable  on  ar.  individual  cluster  size  basis  at  any 
time,  calculations  were  limited  to  size  g=10  after  observing 
virtually  no  numerical  contribution  to  the  results  (i.e.  p,  T, 
X)  above  a cluster  size  of  g=0. 

Initial  conditions  for  both  the  HSC  and  WCK  case  at  273°K 
(saturated  equilibrium  condition)  are  given  in  Table  VI-1.1, 
The  evolution  of  the  cluster  populations  from  the  initial 
condition  {CgJg  Fig.  IV-1.2  is  given  in  Fig.  VI-1.1  through 
VI-1.4.  Recalling  that  the  principal  differences  Ijetween  laoth 
the  classical  IISC  and  the  water  model  lie  in  the  structure, 
the  collision  cross  section  o,  and  the  nearest  neighbor 
parameter  m^^  one  finds  identical  results  for  static  pressure 
and  close  agreement  for  temperature  and  condensate  X.  Also 
the  larger  KCM  cross  section  tends  to  shift  the  hump  in  the 
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TABLE  VI-1.1 

INITIAL  CONDITIONS,  SATURATED  EQUILIBRIUM,  273°K 


Size 

HSC 

WCM 

g 

^g 

% 

1 

C.16000000E+18 

0.16000000E+18 

2 

0.10299099E+15 

0.10299099E+15 

3 

0.37137217E+13 

0.37137317E+13 

4 

0.60239349B+11 

0.14863239E412 

5 

0.47734832E+09 

0.29168346E+10 

6 

0.203534S6E+07 

0.30693876E+08 

7 

0.62537311E+04 

0.23091140E+06 

8 

0.15731362E+02 

0.14066905E+04 

9 

0.34735999E-01 

0.74236071E+01 

10 

0.70293203E-04 

0.35375547E-01 

11 

0.13420118E-06 

0.15648930E-03 

12 

0.24676055E-09 

0.65534307E-06 

13 

0.44375285E-12 

0.26360155E-08 

14 

0.78974221E-15 

0.10297684E-10 

15 

0.14039987E-17 

0.39412301E-13 

distribution  outward  froa  g=4,  HSC  to  g=6  WCM;  the  concentra 
ticR  of  larger  sizes  being  higher  in  the  latter  case.  If, 
for  cxcuaple,  is  used  with  OjjgQ/  the  static  pressure  and 

temperature  calculated  is  about  5%  lower  than  experimental. 
The  disappearance  of  the  smaller  sizes,  akin  to  Ostwald 
ripening  is  also  accentuated  in  the  WCM  case.  Once  the  zone 
of  observed  condensation  (here  x = 7.0  cm)  is  reached  the 
smaller  sizes  (g=2,3,  HSC)  begin  disappearing,  i.e.  dfg<0; 
for  the  WCM  this  applies  to  sizes  through  g=5. 


CLUSTER  SlZEgvWCM 

Figure  VI -1. 3 

variation  of  Cluster  Concentration 
with  Size,  KCM 


2.  Gas  Phase  Metasteible  Equilibrixm 


In  Figure  IV-1.1  the  general  expected  path  of  the  adia- 
batic expansion  up  to  onset  was  described.  Now,  having 
satisfied  all  limiting  constraints  (l-i6)  a final  check  on 
the  solution  is  to  observe  whether  the  adictbatic  path  up  to 
onset  is  calculated  and  mass  balance  maintained. 

Figure  VI-2.1  gives  the  path  calculated  for  both  the 
HSC  and  WCM  with  calculated  results  falling  on  the  adiabat 
for  the  carrier  gas  (Cp/C^  = 1.4);  the  flow  Mach  no.  is  also 
plotted.  The  vapor  pressure  given  is  calculated  from 
Pv  = g^o’^T^g 


Although  the  calculated  pressure  for  4.5<x<6.7  indicates 
that  equilibrium  pressures  are  attained,  the  expanding  mixture 
is  not  in  true  thermodynamic  equilibrium  since  the  gradients 
indicated,  for  example  in  Fig.  VI-1.1,  preclude  equilibrium. 
V7ere  sufficient  time  available  for  the  system  to  evolve  to 
equilibrium  for  4.5<x<6.7,  metastable  cluster  distribution 
as  in  Fig.  IV-1.2  would  exist  (or  conversely,  if  the  rates 
were  sufficiently  rapid) . 

The  results  of  this  calculation  indicate  that  metastable 

equilibrium  is  not  attained  in  the  course  of  the  moist  air 

(8) 

expansion  considered  here  (exp.  405) ' in  spite  of  the 
presence  of  equilibrium  vapor  pressure.  The  vapor  pressure 
up  to  onset  is  a result  of  the  monomer  partial  pressure  which 
follows  the  expansion  hardly  affected  by  cluster  growth  until 
onset  as  is  evident  in  Figs.  VI-1. 2 and  VI-1.4.  Fig,  VI-2.1 


Figure  VI-2.1 

Calculated  Vapor  Pressure  During  Expansion 


indicates  the  vapo;  pressure  from  experimental  inlet  condicion 
to  X = 8.0  cm.  '-diile  interpreting  the  results  after  onset 
(x=6.7),  it  should  be  noted  that  the  temperature  increases  due 
to  the  heat  of  reaction  by  about  24 “K  over  the  dry  isentropic 
temperature.  Over  the  course  of  the  entire  calculation  the 
mixing  ratio  or  specific  humidity  remains  constant  at  Xj^  = 

4.864  g/kg  indicating  conservation  of  mass  throughout.  The 
path  beyond  x=8.0  is  then  again  dominated  by  the  continuing 
expansion  of  the  carrier  ga<’  along  another  adiabat  (GFE,  Fig. 
IV-l.i;  the  process  could  be  interpreted  as  an  expansion  with 
modified  initial  conditions.  That  other  noszle  expansions 

behave  in  thi#  manner  is  evident  by  the  pressure  curve  envelope 
7 8 

after  onset  ' . 

It  has  generally  been  assumed  in  the  past  that  onset  of 
condensation  is  associated  with  collapse  of  metastabie  equili- 
brium; this  would  mean  transition  from  a itetastable  super- 
saturated cluster  distribution.  Fig  IV-1.2  to  a saturated 
cluster  distribution,  III-6.1.  The  condensate  would  then  be 
the  differance  in  the  quantity  g5i9^g  with  the  latter 
distribution  falling  below  the  former. 

Were  Use  reactions  in  our  expansion  sufficiently  fast, 
the  metastable  supersaturated  distributions  of  Fig.  IV-1.2 
would  be  obtained.  It  is  seen,  however,  in  the  non-equili- 
brium evolution  in  Figures  VT-1  through  4 that  the  meta- 
stable  cistribution  of  Fig.  IV"1.2  are  not  attained  - - - 
but  since  the  calculated  vapcr  pressure  equals  that  predicted 
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by  the  equilibrium  isentrope  up  to  onset  point  in  Fig,  VI-2.1 
the  vapor  would  be  assumed  to  be  in  metastable  equilibrium. 

We  see  therefore  that  the  metastable  concept  applies  to  vapor 
pressure  and  supersaturation  but  not  to  cluster  distribution. 

By  the  results  of  Figs.  VI-1. 2 and  VI-1.4  we  can  inter- 
pret onset  of  condensation  either  in  terms  of  monomer  and 
associate  departure  from  metastable  vapor  pressure  with 
disappearance  of  the  monomer  or  alternatively  in  terms  of 
the  entire  cluscer  population  upwards  from  the  dimer.  In 
terms  oj  the  monomer-onset  is  the  'collapse'  of  the  metastcible 
vapcr  pressure,  in  terms  of  the  rest  of  the  cluster  population-- 
onset  is  not  collapse  of  the  metastable  cluster  distribution 
{it  never  having  been  attained)  but  sudden  growth  caused  by 
decrease  in  the  cluster  decay  rates.  This  decrease  is  a 
direct  result  of  the  decrease  in  internal  circulation  frequency 

V. 

It  appears  from  these  calculations  that  the  existance 
of  the  gas  phase  metastable  vapor  pressure,  i.e.  when  J>0, 
is  directly  related  to  both  the  rates  of  growth  and  decay  R^, 
j r"  on  one  hand,  and  to  v on  the  other.  It  also  appears  that 

V may  be  related  to  the  relative  abundance  of  monomer  by 
Equation  v-2,16.  When  sufficient  time  has  elapsed  in  an 
I expansion  for  tho  cluster  concentrations  to  have  grown  due 

I 

j ' to  the  continual  'quench*  and  consequently  the  increased  rates 

iR  , R have  begun  to  affect  the  monomer  concentration,  then 
the  decrease  in  monomer  concentration  causes  the  frequency  v 
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to  decrease  which  further  accelerates  the  process  as  R~  is 
proportional  to  v and  R «R^.  The  net  effect  is  the  release 
of  an  ensrstous  amount  of  energy  in  a short  time  with  the 
temperature  gradient  of  the  flow  changing  to  about  5 x 10^ 
*R/sec.  from  about  -5  x 10^  “K/sec. 

3.  Possible  Water  Cluster  Structures 

The  very  nature  of  the  cluster  theory  originally 
developed  for  moncmolecular  vapors  allows  little  detail  or 
description  of  molecular  structure.  This  atterqpt:  in 
applying  th--  t leory  to  a more  complicated  species  like  H2O 
is  encouraging  especially  since  only  an  average  or  effective 
value  of  pair  interaction  energy  "‘11  is  required  but  it 
raises  the  question  of  the  behavior  of  the  circulation 
frequency  v.  Assuming  that  the  v question  is  resolved  then 
it  appears  that  the  limiting  constraints  developed  here  allow 
extension  of  the  theory  to  other  than  moncmolecular  vapors. 

A composite  summary  of  both  the  haurd  sphere  collision 
and  tlie  water  clathrate  model  is  given  in  Tabic  VI-3.1  and 
Figs.  VI-3.1  and  VI-3.2.  It  should  be  recalled  that  = 1 
for  g = 20  for  both  HSC  and  HCM  which  implies  that  all  mole- 
cules in  this  size  range  of  clusters  are  surface  sites. 

Inasmuch  as  the  theory is  based  on  classical  inter- 
pretation of  xntermolecular  forces,  the  HSC  case  is  in  keeping 
with  this  concept.  A conceptual  view  of  3SC  clusters,  up  to 
g-7  is  given  in  Pig.  VI-3.1  with  an  attempt  to  represent  n^, 
the  number  of  bonds  and  the  number  of  nearest  neighbors 
(taken  here  as  the  average}  . 


! > 


SO 


g = 2. X=  1 
’'2*' 

fX\, 
\ / 

9=3,  X=2 

^3=2 

\ ? 

9 = 4,  X = 2 
>J,=3 

\ > 

9=5,X  = 2*2mj^=  2.02 
175=  4 

The  HCM  requires  some  iiaagination  as  the  linear 


variation  of  is  still  allowed  but  i)  is  l>ased  on  the 

g g 

structural  niodel  of  Fig.  VI-3.2.  Recall  that  as  g-»“(g„); 

here  = 3.9  at  g = 1S5.  Although,  as  pointed  out  earlier 
in  III-4,  numerous  assumed  variations  of  X were  allowed, 

g 

none  other  than  the  present  linear  case  of  X = 2 at  g = 3 to 
X^  satisfied  limiting  constraint  (1-6) . In  the  absence  of 
theoretical  prediction  of  X,  little  more  can  be  said. 

Although  Hg,  WCM  differs  from  n^,  HSC  only  by  approxi- 
mately unity,  this  difference  is  significant  as  appreciable 
energy  is  associated  with  the  bonding  (or  potential  energy) 
term  in  the  reaction  described  by  equation  V-1.11. 

It  is  gratifying  to  find  that  the  parameters  determined 
by  the  series  of  iterations  required  to  satisfy  limiting 
constraints  (1-6)  do,  in  fact,  resemble  those  of  water, 
especially  since  the  value  of  pair  interaction  energy 
HSC  and  KCM,  agrees  with  recent  results  for  an  average  value 
of  2.9  )ccal/mol-bond.  This  average  result  is  deduced  from 
molecular  dynamics  calculations  where  the  detailed  structure 
of  H2O  was  realistically  defined.  The  point  is  also  made  in 
the  molecular  dynamics  study^^®^  that  it  is  legitimate,  to  use 
the  concept  of  a pairwise  additive  potential  in  lieu  of  an 
'exact'  potential  for  water,  prooided  the  pair  functions  are 
effective  pair  potentials. 


J 

\ 
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VII.  COMPARISON  WITH  CLASSICAL  THEORY 


1.  Nucleation  Rates  with  Present  Surface  Tansion 
(3 ) 

Earlier  attempts'  ' to  con^iare  the  predictions  of 
classical  nucleation  theory  with  rapid  nozzle  experimental 

results  were  partially  successful  with  later  studies  reviewed 
(9) 

more  recently.  These  studies  have  in  common  frrst,  the 
application  of  isothermal  equilibrium  nucleation  theory  for 
'birth*  rates  of  critical  size  nuclei  in  highly  dynamic 
supersonic  expansions  (aT/At=-10^®Ksec”^,  or  more) . Second, 
one  is  always  faced  with  assigning  the  surface  tensions  of 
clusters  of  extremely  small  size  of  undetermined  state 
(liquid,  solid).  And,  thirdly,  one  must  choose  a growth  law 
for  clusters  after  their  'birth'  together  with  an  attendant 
condensation  or  'stic):ing'  coefficient. 

(O  Q) 

As  is  evident'  ' ' some  degree  of  success  is  possible 
using  the  classical  approach,  however,  little  progress  can 
be  made  in  interpreting  the  dynamical  aspects  of  phase  change, 
in  understanding  the  metastable  state  or,  in  fact,  predicting 
a priori  the  outcome  of  a new  physical  situation. 

Keeping  in  mind  the  earlier  arguments  of  Sections  V-2 

and  VI-2  with  regard  to  attainment  of  metastable  equilibrium, 

(8) 

the  earlier  calculations'  ' are  repeated  here  using  the  surface 
tension  corresponding  to  the  minimum  in  (Cg)  found  in  Fig. 
rv-1.2.  Plotting  that  value  of  g = g*  at  the  minimum  in 
(Cg)  and  the  corresponding  surface  tension  r(g*),  gives  the 
result  in  Fig.  vn-1.1  with 


Figure  VII-1.1 

Surface  Tension  of  Clusters  if 
Metastable  Equilibrium  were  Attained 


rg*=64.2  + 1.143 (T-200) (dynes/cffi) 


VII-1.1 


An  approximate  relation  for  g*  derivable  from  the  theory 
is 

g*  = 3 - (X„  - X3  - i)  VII-1.2 

which  is  plotted  in  Fig,  VII-1.1.  The  approximation  diverges 

from  exact  values  with  increasing  supersaturation  or  decreasing 

temperature.  At  saturation,  however,  the  parameter  J = In  c^/ 

c,  =0  and  the  approximation  converges  to  the  value  required 
xs 

fay  limiting  constraint  (5) , **  this  is  true  for  both  the 

use  and  HCM! 

Again,  recall  that  our  present  calculations  indicate 
metastable  cluster  distributions  are  not  attained  during  rapid 
expansion  yet  classical  nucleation  theories  predict  the  rates 
for  raetastable  equilibrium  corresponding  to  the  minimum  at 
Cg*.  Thus  it  is  not  surprising  that  repeating  the  earlier 
calculations with  the  value  of  surface  tension  given  by 
equation  VIl-l.l  gives  little  if  any  improvement  over 
previous  results. 

In  Fig.  VII-1.2  typical  results  are  given  for  predictions 

(7  3) 

of  onset  of  earlier  experiments. ' ' ' It  is  possible  to 
iterate  and  converge  numerically  on  experimental  onset  by 
parametric  variation  for,  say,  experiment  405.  Using  the 
iterated  parameters  in  predicting  the  location  of  onset  for 
experiments  416  and  295  with  considerably  different  super- 
saturation  yields  disappointing  results. 
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assSBSSi 


405  BOB,  S,  a =1.0- 


NOTATION 

VOIMER'V 

BECKER,  DORING,  BARNARD>BDB 
VAPOR  PRESSURE  OF  IIQUID«1 
VAPOR  PRESSURE  OF  SOLID  «S 

OBSERVED  ONSET,  EXP  295 — i 


—OBSERVED  ONSET,  EXP  405 
—405  BDB,  S , a * 0.1 
—405  BDB,L,  aO.O 


405  V,  L,  a»1.0 


405  V,  S,  a»l.O 


OBSERVED  ONSET, 
EXP  416 


The  conclusion  one  can  draw  is  that  classical  nucleatiop 


theory  should  not  be  invoked  unless  the  time  scale  of  the 
experimental  situation  is  long  with  respect  to  time  of  evolu- 
tion of  the  metastable  state.  Further  investigation  of  relax- 
ation aspect  of  the  present  theory may  yield  criteria  on 
the  time  of  evolution  of  the  metastable  state. 

2.  Future  Work 

In  order  to  apply  the  experience  gained  with  water  vapor 
to  predicting  other  water  vapor  experiments,  further  investi- 
gation of  the  circulation  frequency  v is  required.  If,  in 
fact,  V = then  other  experiments  should  be  predictable 

from  the  parameters  given  here. 

Reflection  on  the  results  presented  raises  a question 
needing  further  investigation  of  the  net  effect  on  the  heat 
of  reaction  of  the  distinct  approach  used  in  Ref.  (1)  in 
arriving  at  the  dimerization  rate  in  contrast  to  that  for 
larger  sizes.  The  partition  function  is  used  in  the  first 
case,  consequently  accounting  for  the  zero  point  energy  while 
in  the  latter  case  it  is  not  included.  If  the  enthalpy  change 
of  reaction  were  calculated  it  would  be  comptcated  by  this 
fact. 

An  obvious  extension  of  this  work  should  be  to  nozzle 
expansions  of  a monomolecular  gas  such  as  argon  as  well  as 
to  other  species  with  molecular  properties  susceptible  to 
forming  ciathrate  structures,  i.e.  HjS,  H2O,  COj,  SOj,  etc. 
although  argon  itself  might  form  clathrates. 
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A 


*Eq.  8.8+1,  Ref  1 
2nd  term 


ALPHA 

“g 

Surface  to  volume 
parameter 

ALPHA  6 

M 

‘Surface  to  volume 
parameter 

ALPPl 

“g+1 

AREA 

A 

Nozzle  cross  sec- 
tional area 

AREA-MG 

A 

Sub-interval  area 

BIGG 

G 

*Eq.8.9,  Ref  1, 
also  Appendix  4 

BIGM 

M 

*Eq.  8.4,  Ref  1 

BIG«2 

m 

ft 

BK 

k 

‘Boltzman  constant 

BRCFG 

I(3g  - (SgHg)>fg 

C 

s 

*Equil.  concentration 

CF 

•Conversion  Factor 

CHIG 

Vibrational  degrees 

g 

of  freedom 

CHI 

Xg 

* • 

CHIGHl 

Xg_l 

CHIGPl 

Xg+1 

CHIl 

Xl 

•Syinbols  for  equilibrium  solution  that  differ  from  nonequili- 
brium solution. 


CKP 

•Pressure  conversion 
factor 

COND 

X 

Condensate  mass 
fraction  in  g/kg 

CP0 

s 

Specific  heat 

CUS 

•Under-saturated 

9 

equilibrium  concen- 
tration 

CV0 

Cv 

Specific  heat 

ClINF 

Saturated  equil.  C^ 

C10 

C. (frozen  equil., 
see  Appendix  2) 

Cl 

‘^l 

•Equil.  monomer  cone. 

C10 

C-^  (initial) 

•Initial  monomer  cone 

Cll 

<=11 

•Eq.  8.2,  Ref  (1) 

C12 

<=12 

•Eq.  8.4,  Ref  (1) 

CIG 

•Eq.  8.4,  Ref  (1) 

CIKE 

•Eq.  III-8.7,  2nd 
term 

C2KE 

•Eq.  III-8.7,  3rd 
terra 

CGKE 

•Eq.  III-8.7,  4th 
terra 

CGVIB 

•Eq.  III-8.7,  5th 
terra 

CGHK 

•Eq.  III-8.8,  6th 
terra 

DAOX 

dA/dx 

Slope,  nozzle  area 

DAOXMG 

Sub-interval  dA/dx 

DAHG 

Sub-interval  dA 

DDAOX 


d (Oadx) 


I 


lr 


DUU 

DELEP 

DELEP0 

DELTAP 

DELTAT 

DELTAX 

DELZI 

DELZET 

DENOM 

DENOHl 

DENOM2 

DENIS 

DFG 

DFGA 

DFGEX 

DFGHl 

DFGPl 

DF0 

DF10D 

DFG12 

DFGIEX 

DFG2 


d^u 

See  Eq.  V-1.4 

Ac” 

g 

Eq.  V-1.2 

°static~^isen tropic 

^static  — T.  . 

reaction  ^se»tropic 

Experimental  data 
interval 

AC 

S - Vi 

m 

* • 

Denominator/  Eq. 
V-2.7 

First  term,  DENOM 

Second  term,  DENOM 

DENOM,  no  reaction 
(isentropic  case) 

Eq,  V-1.13  (total) 

<“g 

Eq.  V-1.12  (reaction) 

Eq.  V-2.11  (expansion) 

Eq.  V-1.13 

^g^i 

Eq.  V-1.13 

Eq.  V-2.12 

Initial 

Initial  value,  no 
reaction 

dfj 

Eq.  V-2.15 

^1 

Eq.  V-2.11 

^2 

Eq.  V-2.13,  g=2 
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DFG2M1 

DFG2,  previous  time 
increment 

DLDG 

“x 

Eq.  III-4.1 

DPISEN 

dPi 

Eq.  V-2.10,  no 
reaction 

DPSTAT 

dp 

Eq.  V-2.10  or  V-2.9 

DQ 

dq 

(dq/dt) dt 

DQDT 

dq/dt 

Eq  VIII-3 

DQDX 

dq/dx 

(dq/dt) /u 

DTISEN 

dt. 

Eq.  V-2.7,  no  reaction 
term 

DTM 

dt 

in  seconds 

DTSTAT 

dT 

Eq.  V-2.? 

DU 

du 

Experimental  velocity 
increment 

DDDXHG 

du,  calculating 

DXCH 

dx 

Increment  of  length 
along  nozzle  eucis 

D1 

Eq.  VIII-3 

D2 

M 

D3 

n 

EPSIF 

^Iteration  limit, 
statement  240 

STAG 

\ 

Eq.  III-8.2 

FGIM 

f- 

Hon-stationary 

9 

concentration  (cm”'’) 

FI 

•Eq.  (8.8)+l,  Ref  (1), 
3rd  term 

FIRST 

Eq.  VI-2.1,  1st  term 

FIRST 

•Eq.  III-8.7,  1st  term 

F0TM 

f 

o 

C2urrier  gasconcen- 
tration  (cm”"’) 

F0IK0 

initial 

FTHETA 

F(e> 

‘From  Eq.  8.6,  Ref 
(1) 

FRT 

F(6) 

*Eq.  8.15,  Ref  (1) 

FTHETG 

F(g,0) 

*£q.  8.6,  Ref  (1) 

F10D 

f , , no  reaction, 
iiitial 

FITM 

Non-stationary 
monomer  concentration 

F1TK0 

fj^,  initial 

FICIR 

f same  temp. 

Eq.  VI-2.1 

FITERM 

Eq.  V-2,7,  2nd  term 

F2TM 

^2 

Non-stationary  dimer 
concentration  (cm“3) 

GAMA 

Y 

Numerical  constant, 
see  App  3 

GAMAS 

r 

Surface  tension, 
dyne/cm 

GAHAO 

dimer 

GEAR 

i 

Average  cluster  size 

GGPKG 

X 

Eq.  V-2.12 

GUlOA 

X 

g 

Coordination  number 

GHAX 

* = MAXG,  below 

GINF 

9- 

numerical  infinity 
for  cluster  size 

GRATIO 


G<e)/G(g,3) 


>L  /L  , pg  1282,  Ref 


GTHETA 

G{g,6) 

*Eq.  8.9,  Raf  1 

LG 

g-1 

H 

g+1 

HAXG 

tipper  limit  on  actual 
calculation 

MAXGPl 

MAXG  plus  1 

MG 

Siibinterval  index 

NG 

Calculating  limit 

NGAHA 

Experimental  interval 
divisor 

NSTEP 

Stepping  index 

K0 

Index  locating  start 
of  calculation 

PAIR 

Initial  partial 
of  an  Isentropic 
pressure,  cm  Hg 

PCOEF 

Conversion  factor 

PI 

3.14X59 

PINF 

Ps 

Equilibrium  saturation 
vapor  pressure 

PISEN 

Pi 

Isentropic  static 
pressure  (no  reaction) 

PPAR 

Pv 

Partial  pressure  of 
vapor 

PPAR0 

(p^)=  initial 

PRES 

*Calc.  pressure,  cm. 
Eg. 

PSTAT 

p 

Static  pressure  of 
ceurrier  gas  vapor 

PVAP 

Partial  pressure  at 

initiation  of  program, 
ref.  only 
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I 


! 

il 


1 

? 


1 


P(»l 

static  pressure 
at  inlet  conditions 

QOFX 

^(x) 

Z^dq,  See  DQ 

RDTM 

dt 

Time  increment 
(corrected) 

REF 

Eq.  V-2.15 

RMATG2 

dfg/dt 

Eq.  V-1.12 

RMAT0 

Eq.  V-2.12 

RMAT22 

Eq.  V-1.12,  g=^2 

EMATIO 

Eq.  V-2.H,  F, 
isentropic 

EMX'rj 

% 

Table  V-'l.l 

RMIUFl 

Table  V-1.1 

RMI2BI 

R- 

Table  V-1.1 

RMI3U 

Table  V-1.1 

RMI4U 

‘‘4 

Table  V-1.1 

RPLUS 

Table  V-1.1 

RPLGPl 

Table  V-1.1 

RTAR 

Eq.  V-2.5  (see 
conraent  in  prograua) 

RTUU 

Eq.  V-2.5  (See 
comment  in  prograua) 

R0 

Mean  molecular  radius 
carrier  gas 

R1 

r. 

Meeui  molecular  radius 

X 

vapor 

S 

s 

Saturation  ratio 

SIGMAl 

^i,g 

Eq.  II1-1.5,  collision 
cross  section 
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SI601 

“0,1 

Collialon  cross- 
section,  carrier 
gas  - Btoncaier 

SZG02 

“0,2 

Collision  cross- 
section,  C2trrier 
gas  - dimer 

SIGH 

Collision  cross- 
section,  mononer- 
monomer 

SIG12 

“1,2 

Collision  cross- 
section,  monomer 
dimer 

SIGIG 

“i,g 

*See  SIGNAL 

SIGIGI 

“i,  g-i 

*See  SIGHAl 

SIG101 

Collision  cross 

r 


SIJIADX 

dWdx'‘ 

Slope  of  DAOX 

SLDU 

d^u 

SLU 

du 

Slim 

du  at  x-1 

SMGDFG 

Egdfg 

Alternate  expression 
for  Eq.  V-1.5 

SMNGRG 

Eg.  V-1.6,  1st  term 

SRDTM 

t 

Z dt 

SUH 

log  V 

Eg.  VI-2.1 

SOMCHI 

See  Eg.  1II-8.7 

SUKCG 

EC, 

« II  Tl 

SUHGCG 

Egc, 

Eq.  8.17,  Ref  {!) 

SDHDAX 

Zi (OAOX) 

See  top  3,  multipoint 
derivatives 

5.^ 


SUMDDU 

Zti  (DD) 

SUMDU 

2DU 

« 

SUHBRC 

SUMGFG 

“c, 

• 

SUMGG 

SUMLCG 

SUMLNW 

S{(uywj 

SXER 

d/dx  (XER) 

TDK 

T 

TDKIS 

’’i 

TDK0 

TDK00 

'^OO 

THETA 

e 

U 

u 

UFAC 

UHG 

HG 

“g 

WG3 

“3 

HINF 

vsm 

u® 

"oo 

See  App  3,  multi- 
point derivatives 

n 

Eq.  V-1.8,  Z{  } 


See  1st  term,  Eq. 
III-8.7 


See  eq.  8.13,  Ref  (1) 
Appendix  VIII-4 
Temperature  ®K 


Temperature, 
isentropic  ®K 

Temperature  at 
initiation  of 
calculations 

Temperature  at 
inlet  conditions 

6 = u'j^j^/JcT 

Flow  velocity 
u®ii  = UFAC  X 


Flow  velocity,  sub- 
interval 


Eq.  8.3,  Ref  1 
* for  g=3 

" for  g=» 


Pair  interaction 
energy  carrier  gas- 
carrier  gas 
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Pair  interaction 
energy  carrier 
gas  - monomer 

Pair  interaction 
energy  monomer  - 
monomer 

distance  on  nozzle 
centerline 

multiplier,  see 
program 

Initial  XER 

Coq>uter  index 
(note  X3) 

*Eq.  8.5,  Ref  (1) 


*See  eg.  S.IO,  Ref 
(1) 

•Bq.  8.9,  Ref  (1) 


XLGTDK 


XIOSCG 


Eg.  y-2.14 

*Eg.  III-8.7 

*See  eg.  8.10,  Ref 
(1) 

*See  eg.  8.10,  Ref 
(1) 


Eg.  V-2.14 
Sg.  IV-1.2 


XLRIIO 


*Pg.  1282,  Ref  (1) 


J,-4J 


Specific  humidity 

Molecular  mass, 
carrier  gas 

Molecular  mass, 
vapor 

Reduced  mass, 
carriar  gas  - cluster 

Reduced  mass, 
monomer  - cluster 


Reduced  mass, 
monomer  - dimer 

Reduced  mass, 
carrier  gas  - dimer 


NONDAT 


Input  data  file 


Input  data  file 

Circulation  frequency 

Eq.  V~2.7,  numerator 

Eq.  V-2.7,  1st  & 2nd 
terms  of  numerator 

Eq.  V-2.7,  3rd  term 
of  numerator 

Eq.  V-2.7,  no 
reaction 

Output  file 

•Eq.  8.13,  Ref  (1) 

Monomer  rotational 
mode 

Eq.  V-2.11, 
collision  number 
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FORTRAN  PROGRAM  FOR  CLUSTER  DISTRIBUTIONS  AT  SATURATED  EQUILIBRIUM 

The  solution  given  is  for  water  vapor  in  air  but  it  can  be 
applied  to  any  mixture  of  a vapor  and  a carrier  gas.  Inputs 
required  are;  T.  g„,  e,  r^,  r^^.  u°j^j^,  m^,  p^, 

Cj^,  and  the  nomenclature  used  being  given  in  Appendix  1. 

The  iteration  involved  is  given  in  Table  III-2.1  and  Flow  Chart 
Vlll-2.1. 

A term  in  the  program  but  not  yet  described  is  UFAC.  This 
is  a multiplier  which  facilitates  the  alteration  of  u®  faTOID 
the  value  of  = ISOOK  (erg/®K)  (or  2.98  K cal/mole  assumed 

initially). 

A further  remark  concerns  According  to  the  cluster 

theory hard  sphere  collision  diameters  are  to  be  used. 

The  value  used  here,  however,  is  not  = 2.00r-  but  rather 
Oil  = 2.05rj^.  This  value  was  obtained  by  cross  iteration 
between  the  present  equilbrium  solution  and  the  non-equilibrium 
solution  of  Appendix  3.  The  value  of  o^^  was  adjusted  until  the 
cluster  reaction  rates  beceune  stationary  to  about  5 digits  when 
the  equilibrium  distribution  of  cluster  sizes  was  fed  into  the 
non-equilibriiOT  program.  Vapor  pressure  data^^*  for  liquid 
water  when  fitted  by  least  squares  approximation  yields  p^  = 
9,3286  - 2367. /T  while  the  expression  for  adiabatic  expansion 
of  moist  air  from  an  initial  temperature  of  295®K  and  specific 
humidity  of  about  5g/kg  is  P «•  0.425  (T/273)^'®.  The  pressure 
conversion  factor  cm  Hg  to  dynes/cm  is  1.33  x 10  . 


Once  ini.tiai  conditions  have  been  selected,  the  iteration 
in  Table  IH-2.1  taay  proceed  with  the  objective  of  satisfying 
limiting  constraints  (1)  through  (5) ; condition  (6)  is  dealt 
with  later  in  Appendix  3.  For  each  cluster  size  from  g = 2 
upwards  to  an  eurbitrary  upper  limit  equation  III-5.1 

must  be  solved  for  the  concentation  (cm“^) . Also  calculated 
is  The  approach  of  to  zero  as  g approaches  numerical 

infinity  can  be  checked  numerically  in  file  F2  or  plotted 
separately  for  inspection. 

To  calculate  a super  or  undersaturated  cluster  distribution 
one  must  first  calculate  the  saturated  distribution  at  the  seune 
teaq>erature  to  obtain  the  value  of  Cj^  in  statement  915  for  CUNT. 
Statonent  800,  P “ must  be  deleted  and  statement  1010  deleted. 
All  other  parameters  remain  the  same  as  for  the  saturated  solu- 
tion at  that  temperature. 
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APPENDIX  3 


PORTRAN  PROGRAM  FOR  NON-EQUILIBRIUM  CLUSTER  DISTRIBUTION  DURING 
AN  ADIABATIC  EXPANSION 

The  calculation  given  here  is  for  the  expansion  of  water 

vapor  in  air  in  a supersonic  nozzle  of  known  geoaetric  profile. 

11  8) 

Solution  of  the  one  dimensional  equations  of  motion  ' ' ' yields 
the  flow  velocity  u which  is  required  as  an  input  datum.  The 
same  solution  gives  the  cross  sectional  area  A and  its  rate 
of  change  along  the  flow  axis  dA/dx.  Also  required  are  the 
initial  equilibrium  cluster  distribution  from  Appendix  2 and 
the  cluster  properties  such  as  radius  r,  pair  interaction 
energy  . 

The  computational  scheme  or  flow  chart  is  given  in  Figure 
VIII-3.1  and  the  relation  of  the  iteration  time  interval  to 
data  increments  given  in  Figure  VIII-3.2.  The  method  is  to 
divide  the  data  interval  (6x  = 1 mm)  into  sub  intervals  (Ax/t) 
small  enough  to  allow  a non-oscillating  solution  by  Eulerian 
numerical  integration: 

f (x  + Ax/y)=  f(x)  + d f/dx(Ax/y)  VlII-3.1 

However,  the  fewer  the  subdivisions  the  shorter  the  computation 
time.  The  smallest  value  of  y giving  stability  was  found  to 
be  about  5000  and  this  value  already  brings  the  calculation 
interval  down  to  the  order  of  the  molecular  collision  time. 

An  extensive  effort  was  necessary  to  find  a method  to  reduce 
the  eunount  of  computer  time  involved.  Additional  complications 
in  that  the  collision  time  increases  during  the  expansion  process 


and  the  numerical  interval  in  Equation  VlII-3.1  must  be 
significantly  longer  than  the  collision  intcirval  to  be  con- 
sistent with  equilibration  of  translational  energies.  An 
approximation  technique  yielding  stable  numerical  results 
was  developed  after  many  trials.  After  dividing  a data  interval: 
AX'iX/y.  with  a value  of  y around  5000,  a multiplier,  XER  for 
example,  XER=2.0  is  selected  arbitrarily.  This  multiplier 
has  the  effect  of  lengthening  the  calculation  interval  because 
the  second  term  in  Equation  VIII-3.1  becomes 

Ax/y  (XER)  d F/dx  VIIl-3.2 

1 

To  reduce  the  number  of  calculation  intervals  within  the  { 

original  data  interval  Ax,  an  end  condition  ESTOP  is  introduced 
where  NSTOP  = Y/XER.  If  XER  is  2.0,  the  derivative  d F/dx  is 
thus  applied  over  a length  for  each  y interval,  but  the 
calculation  for  a given  data  interval  is  terminated  after  y/2 
sub  intervals.  To  take  account  of  the  increase  in  collision 
time  due  to  expansion,  the  collision  time  for  the  g = 10 
cluster  is  calculated  at  the  beginning  and  end  point  of  the 
overall  calculation  (an  increase  from  about  .5  nsec  to  4nsec) 
and  a linear  fit  to  XER  is  determined.  An  increment  SXER  is 
added  to  XER  and  NSTOP  recalculated  before  calculating  over 
the  next  data  interval.  Thus,  a crude  correspondence  between 
the  collision  interval  and  the  calculation  interval  is  main- 
tained and  the  amount  of  computer  time  is  reduced  since  NSTOP 
is  decreasing  continually. 

In  as  much  as  the  calculation  proceeds  at  near  the  molecular 
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collision  interval,  interpolation  of  the  experimental  input 
data  and  the  calculated  flow  velocity  u is  required.  Multi- 
point (seven)  derivatives  are  taken  over  the  original  numerical 
data  and  velocity  then  applied  each  V subinterval.  Thus,  the 
original  data  for  nozzle  slope,  DAOX(J)  (INDEX  J increments 
each  m.ffl.)  becomes  DAOXMG  for  each  y interval.  The  expansion 
terms.  Equation  V-2.11  and  V-2.12  are  then  corrected  each  y 
interval. 

Additionally,  after  calculating  all  rates  up  through  g^^^ 

at  any  subinterval,  the  rates  for  cluster  size  (g_,  +1)  must 

be  calculated  to  include  the  influence  of  on  fg.  In  this 

calculation  the  influence  of  is  ignored  so  that 

’max  + 2 


df„ 


’max+l 


d (£n  V) 


gtl  - «gtl>  - fgtl  3t 
rather  than  Equation  V-2.13.  This  result  is  correct  within 
a few  percent  as  long  as  ^o+2^^^g+l'  calculation  appears 

after  statement  525  in  the  program. 

An  exact  calculation  over  all  y subintervals  and  data 
intervals  involves  about  16  hours  on  the  PDF  15/40.  With 
X£R=2.0,  the  time  is  reduced  to  about  4 hours  while  on  using 
the  SXER  technique  the  time  is  about  1.5  hours.  This  is  a 
valuable  saving  in  computer  time.  No  significant  difference 
in  the  results  is  obtained  between  these  three  methods. 

The  initial  debugging  and  progreun  development  of  this 
program  and  the  equilibritsn  program  was  done  in  an  interactive 
mode  over  several  years  on  a H-635  and  DTSS  (Deirtmouth  Time 

lie 
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Sharing  System) . The  results  described  would  have  been 
impossible  to  achieve  without  this  computer  service  where  50 
or  more  short  runs  could  be  accomplished  in  one's  own  office 
in  two  hours.  Once  debugged  the  programs  were  shifted  to  the 
POP  15/40  for  the  long  runs.  The  accuracy  could  also  be  checked 
as  7-8  digits  are  significant  on  the  H-635  while  only  5-6  digits 
are  significant  on  the  POP  15/40.  Using  double  precision  on 
certain  quantities,  closely  similar  results  were  obtained  with 
the  two  computers. 

Another  feature  that  should  be  mentioned  is  the  restart 
capability.  If  for  some  reason  a run  is  terminated  before 
completion,  it  can  be  restarted  with  loss  of,  at  most,  less 
than  one  data  inter /al  of  computation.  The  restart  is  not 
precise  since  only  average  values  of  the  nozzle  geometry  are 
used  (see  RESTART  PATCH) , however,  if  for  some  reason  long 
runs  become  infeasible,  the  restart  information  could  be  fed 
into  the  initial  file  (OUTFIL  SRC)  each  time  NSTOP  is  reached 
(Statement  514  + 1) . 

When  NSTOP  is  reached  the  value  of  all  cluster  reaction 
terms  (R*,  r”,  “^^g'  flow  properties 

(p/Poi»  etc.)  are  written  out.  Isentropic  properties 

calculated  by  suppressing  the  reaction  term  of  the  expression 
for  dT.  (Equation  V-2.7)  are  also  reported.  The  partial 
pressure  of  the  vapor  obtained  by  summing  the  partial  pressures 
of  each  cluster  size  is  reported  as  PPAR. 
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CORRIGENDA  TO  REF.  (1) 
TRANS.  FAR.  SOC.  65,  1267,  1969 


Equation  3.4  = 9 “l  ^ 


,2/3 


Equation  7.1+1  (y^)^  (e/3g) 


3g 


»l/2 


Equation  8.5 
Equation  8.8  G (3,0)=2.O (e®  - 1)/(1  + 29) 
Equation  8.9  G(g>3,9)  = {exp[9  (^-i)  l"sxpt9 


(1  + ^9+!^  hh 


Equation  9.1  Rj^  = Z (R  - R^)  + 2 {R2  ~ R^2^ 
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